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The paper is devoted to the analysis of a spectral collocation method for approximating 
stationary Navier-Stokes equations, governing the flow of a viscous incompressible fluid in 
a domain 0 of (R 2 or OR 3 

, (1.1) - v Au + grad p + (u.V)u = f in 0 

div u = 0 in Q , 

I provided with Dirichlet boundary conditions 

(1.2) u = g on 90 

In this system , the data are the density of body forces f, the velocity on the boundary g and 
the kinematic viscosity v > 0 ; the unknowns we want to approximate are the velocity u and 
the pressure p inside the domain 0. The discrete solution is sought in a space of polynomials 
of high degree on 0 , and the equations (1. 1 )( 1 .2) are verified in a finite number of points, 
called collocation nodes. We refer to [VGH] and to [CHQZ] for a detailed bibliography about 
such methods. As far as the theoretical analysis is concerned, we limit ourselves to the test 
domain 0 = ]-1 ,1 [ 2 . In this square, the collocation points form a cartesian grid : their 
coordinates belong to the set of the nodes of a Gauss-Lobatto quadrature formula. However, 
the numerical results prove that there is no difficulty to extend the method to 
three-dimensional and/or curved domains. 

The method we analyse in this paper has the following features : 

1) Only one grid is involved in the algorithm : indeed, all the equations of (1.1) will be 
satisfied in the same nodes. We refer to [MMo] and [BM2] for other collocation methods for 
the Navier-Stokes equations, which are called staggered grid methods. 

2) Due to this unique grid, the velocity and the pressure are approximated by polynomials 
of the same degree. Then, it turns out that there exist some spurious modes for the 
pressure, i.e. some polynomials the gradient of which cancels at the collocation nodes; this 
fact has first been pointed out in [Mo]. In order to obtain the convergence of the pressure, it 

' is necessary to choose a suitable discrete pressure space which does not contain these modes 

but retains good approximation properties. We shall use the space already proposed in 

* [BMM J[Me] or [BCIi]. 

3) The grid can be built from the Gauss-Lobatto quadrature formula associated with any 
family of Jacobi polynomials. Here, we shall treat two special cases. The first one is the 
case of Legendre polynomials; it is the simpler one, since its analysis involves the standard 
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variational formulation of the Navier-Stokes equations in standard Sobolev spaces. 
However, we also consider the case of Chebgshev polynomials; indeed, the nodes are then 
images by the cosine function of equidistant points, so that the use of the Fast Fourier 
Transform allows for a less expensive computation of the derivatives or of the nonlinear 
terms. This last method is numerically cheaper, but its analysis involves a non trivial 
formulation of the equations in appropriate weighted Sobolev spaces with the Chebyshev 
weight, the properties of which are given in [BM 1 ]. 

The method we present has already been studied in a simpler case : for the Stokes 
problem that one obtains by neglecting the nonlinear terms in (1. 1 ) 

(1.3) -vAu + gradp = f in 0 , 

div u = 0 in 0 , 

when it is provided with homogeneous Dir ichlet boundary conditions 

(1.4) u = 0 on 50 

We only recall the results of [BMM] in the Legendre case, of [BCM] in the Chebyshev case. 
Then we extend them to the Stokes problem with non homogeneous boundary conditions, by 
using the lifting of polynomial boundary data of [BM 1 ]. Finally, the method is applied to the 
full Navier-Stokes equations ; the nonlinear terms are handled in a pseudospectral way. The 
justification follows from the discrete implicit function theorem of [BRR], in a slightly 
different form due to [C] (see also [CR ] ) . The error estimates we obtain between the 
discrete solution and the exact one, when it is assumed to be smooth enough, are the same as 
for the Stokes problem. This theoretical justification is completed by numerical 
experiments, achieved in the case of a three-dimensional domain with Chebyshev nodes, and 
a description of the algorithm involved in the code is given. We refer to [Me] for more 
details on the numerical implementation. 

An outline of the paper is as follows. Section II is devoted to the definition of the 
discrete approximation spaces and collocation problems, first in the homogeneous case, then 
in the inhomogeneous one. In Section III, we recall the convergence results of [BMM] and 
[BCM] for the Stokes problem in the homogeneous case, then we complete them for 
inhomogeneous boundary conditions. In Section IV, the analysis is extended to the 
Navier-Stokes equations, in both cases of homogeneous and inhomogeneous boundary 
conditions. Finally, in Section V, the techniques required by the numerical implementation 
of the method are presented and examples of numerical results are given. 
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Notation 

The norm of any Banach space E is denoted by ||.|| E . For any pair (E,F) of Banach 
spaces, £(E,F) represents the space of continuous linear mappings from E into F. We mean 
by A ® B the tensorial product of any sets A and B in a Banach space, while A® 2 is the 
tensorial product of A by itself. In all that follows, c, c' ... are generic constants, 
independent of the discretization parameter. 

In Sections II to IV, we shall work in the square 0 = ]- 1 , 1 [ 2 . Let us precise some 
notation about this domain. The generic point in 0 will be denoted by x = (x,y) (or 
sometimes by (x, , x 2 )). We introduce thecornersaj , J € Z/4Z, of 0 (where a Jt1 follows 
aj counterclockwise), and call Tj the edge with vertices aj_, and aj ; for any edge Tj , J € 
Z/4Z, nj is the unit outward normal to 0 on Tj and tj the unit vector orthogonal to rij , 
directed counterclockwise. 



Sobolev spaces 

For any domain A in IR d and for any real number s > 0, we use the standard hilbertian 
Sobolev spaces H S (A), the norm of which is denoted by ||.|| s A . On the square Q , we shall use 
the scalar product 
(1.5) (<p,<j<) = \ n <p(x) «|»(x) dx . 

We also recall that, for any integer m ^ 1 , the semi-norm 
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0.6) =< 0 lla m «/ex J a u m - J llI. 0 ) m 

is a norm on the space H™ (Q) of the functions of H m (0) which vanish on the boundary 90 
together with all their derivatives up to order m- 1 (the traces being defined in the sense of 
[LM]). The dual space of H™(0) will be denoted by H' m (Q), and it is standard to note that 

(1.7) H _1 (0) = { f + 0g/9x + 5h/3y, (f,g,h) € L 2 (Q) 3 } . 

The space of functions in L 2 (Q) with a null integral is noted Lg(Q). 

Next, we recall some basic material about weighted spaces of Chebyshev type (for 
further details, see e.g. [CHQZ][BM 1 ][M2]). If p(U = ( 1 -C 2 )' ,/2 denotes the Chebyshev 
weight on the interval A =]- 1 ,1 [, let 

L 2 (A) = { (p : A -> 1R ; <p 2 (C) p(C) dC < +oo } 

be the Lebesgue space associated with the measure p(C) dC , provided with the norm 

(1.8) ll«Pll 0 ,.A = (I-,-P 2 (Op(C)dC) ,/2 • 

The weighted Sobolev spaces are defined as follows : for any integer m ^ 0, H™(A) is the 
subspace of L 2 (A) of the functions such that their distributional derivatives of order < m 
belong to L 2 (A) ; it is a Hilbert space for the inner product associated with the norm 

a?) Mu* = ( l:.. i< iA r' 2 . 

where 

(MO) l«f>l k , p . A = ||d k «p/dC k || 0 . p .A • 

For a real number s = m + a, m € IN, 0 < a < 1 , we define H*(A) as the interpolation space 
between H™ +1 (A) and H™(A) of index 1-a (cf. [LM]); we denote its norm by IUI M A . 
Finally, we can apply a rotation and a translation to define similar Sobolev spaces on any 
segment of length 2 in IR 2 . We use the same notation as before to indicate them , as well as 
their norms. 

The Chebyshev weight on the square 0 is defined as cj(x) = p(x) p(y). Let 
L 2 (Q) = { (p : 0 -> IR ; J 0 <p 2 (x) cj(x) dx < +oo } 
be the Lebesgue space associated with the measure w(x) dx, provided with the inner product 
(Ml) (ip,ip) u = j ft <p(x) «p(x) w(x) dx 

and the norm || || 0uft = (. .-)J[ } /2 - Next, weighted Sobolev spaces are defined as follows : for 
any integer m > 0, H^(Q) is the subspace of L 2 (0) of the functions such that their 
distributional derivatives of order < m belong to L 2 (G); it is a Hilbert space for the inner 
product associated with the norm 
(112) II ^ Hm.w.O = ( = o I Ik ,u,0 ^ 


where 
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(M3) l<pl k , u , 0 = ( Lj s0 lla%/a x j ag k - j || 0 2 u<0 ) 1/2 . 

For a real number s = m + cr,m€(N,0<cJ< I, we define H^(Q) as the interpolation space 
between H™ +, (Q) and H™(Q) of index 1 -a; we denote its norm bg ||.|| su 0 . Finally , for any 
integer m > 1 , we consider the closed subspace of the functions of H™(Q) which vanish on 
the boundarg 30 together with all their derivatives up to order m-1 (the traces being 
defined in the sense of [LM]); this space, denoted bg 0 (O), is the closure of J2>(0) for 
the norm of H™(Q) (see [BM 1 , Prop. 11.9]). Due to the Poincare inequalitg, an equivalent 
norm on H™ 0 (Q) is the semi-norm |.| m u0 . The dual space of H™ 0 (ft) will be denoted bg 
H‘ m (0); if the space L 2 (Q) is identified to its dual space, we have for instance 
(I.M) H- 1 (0) = {f + eg/a X + ah/ag,(f,g,h)eL 2 (0) 3 } . 

We also introduce the space L^ 0 (O) of functions q in L 2 (Q) such that J n q(x) u(x) dx is 
equal to 0. 


Variational formulations 

In order to treat the Legendre and Chebgshev approximations simultaneouslg, we 
introduce a letter A which is L in the Legendre case and C in the Chebgshev case, a 
parameter oc equal to 0 in the Legendre case and to —1/2 in the Chebgshev case (this is the 
power of (1-C 2 ) involved in the corresponding weight). For instance, the sgmbol H^(Q) 
stands for the space H s (0) in the Legendre case (A = L, a = 0) and for the space H^(0) in 
the Chebgshev case (A = C, oc = - 1 /2). 

To write appropriate variational formulations of equations (1. 1 ) and (1.3) , we first 
consider the boundary condition (1.2). Let us assume that the function g is such that the gj = 
g lrj , J € Z/45Z, satisfg 

(1.15) g 0 € H^-“ )/2 (rj) 2 , J € Z/4Z , 

(1.16) cz/42 |pj flj • n j = 0 • 

Assume moreover, in the Legendre case, 

(U7) L lo CSj(aj - t tj) -g j+1 (aj + t -C j + ,)] 2 t' 1 dt < +co ,J€Z/4Z , 

and, in the Chebgshev case, 

( 1. 1 7) c gj(aj) =g J+1 (aj) ,J€Z/4Z . 

Then, there exists [G, Thm 1.5.2.3][BCM, Thm II 1.2] a function u b in H^O) 2 satisfging 

(1.18) div u b = 0 in 0 , 

(1.19) Ub = 9j on ^*j .J€Z/4Z • 
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Next, we define the bilinear form a A on H A (Q ) 2 x H A 0 (Q ) 2 by 

( 1 . 2 0 ) L V u € H 1 (0 ) 2 , V V € Hq( 0) 2 , a^u.v) = v (grad u , grad v) , 

( 1 . 2 0 ) c V u € H^(0) 2 , V ve H^ 0 (O) 2 , 

a c (u,v) = v | 0 (grad u)(x) (grad (vw))(x) dx . 
Clearly, for any f in H A '(Q ) 2 and any g satisfying (1.15), (1.16) and (1. 1 7) A , problem 
(1. 3)0.2) is equivalent to the following variational one : Find a pair (u,p) in 

H a (Q ) 2 x L ao (0) , with u-u b in H A 0 (Q ) 2 , such that 
( 1 .2 1 ) A Vv<eH A 0 (Q) 2 , a A (u,v) + (v , grad p) A = (f,v) A , 

V q € L A (Q) , (div u , q) A - 0 . 

In the Legendre case, it is well-known that problem (1.2 1 ) L admits a unique solution. In the 
Chebyshev case, it is also known [BCM, Thm III. 2] (but may be not so well) that problem 
(1.21 ) c admits a unique solution. In both cases, the solution satisfies the stability estimate 

(1.22) INI, A 0 + II Pll 0 ,A.n ^ C ( II^H^Vo ) 2 + ^ J 6 2/42 II 8 jll( 1 -«)/ 2 ,A,rj ) 

Of course, if the solution (u,p) of (1.21 ) L belongs to ( Q ) 2 x l 2 (0), then the pair 
(u , p - ( 1 /n 2 ) J 0 p(x) cj(x) dx) is the solution of problem ( 1.2 1 ) c . 

As far as the Navier-Stokes equations (1. 1 )(l.2) are concerned, for any f in H^Vo ) 2 
and any g satisfying (1.15), (1.16) and (1. 1 7) A , they admit the following variational 
formulation: Find a pair (u,p) in H A (0 ) 2 x L 2 0 (O) , with u-u b in H A 0 (O ) 2 , such 
that 

(1.23) a V v € H ao (Q) 2 , a A (u,v) + (v , grad p) A + ((u.V)u , v ) A = (f,v) A , 

V q € L A (0) , (div u , q) A = 0 . 

In the Legendre case, this problem admits at least one solution. If this solution belongs to 
H,1(0 ) 2 x L 2 (0) , it is also a solution of ( 1 . 23 ) c , up to an additive constant on the pressure. 
More details will be given in Section IY. 
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IL The collocation problems. 

We begin bg introducing the collocation framework, especially the collocation grid. 
Then, we present the collocation discretization of the Stokes and Navier-Stokes equations 
provided with homogeneous boundary conditions. That leads us to define suitable discrete 
spaces of pressures. Finally, we can extend the collocation problems to the case of 
inhomogeneous boundary conditions. 

11.1 ■ The collocation .framework. 

Let us introduce some monodimensional notation. For any nonnegative integer n, 
P n (A) denotes the space of restrictions to A = ] — 1 , 1 [ of polynomials of degree < n. We 
shall use two families of orthogonal polynomials on A : 

1) the Legendre polynomials ( L n ) n?0 , which are orthogonal for the measure dC , 
normalized by the following condition : the Legendre polynomial L n , n > 0, is of degree n 
and satisfies L n (± 1 ) = (± 1 ) n ; 

2) the Chebyshev polynomials ( T n = cos (n Arccos C) ) n?0 , which are orthogonal for the 
measure p(C) dt; ; of course, the Chebyshev polynomial T n , n ^ 0, is also of degree n and 
satisfies T n (+ 1 ) = (± 1 ) n . 

In order to have a unique notation in the Legendre and Chebyshev cases, we introduce, for 
each real number <x > -1 , the Jacobi polynomials ( J“ ) n?0 wh ’ ch are orthogonal for the 
measure ( 1 -C 2 )* dt;. Since J* is of degree n and such that 
J“(±l) = (±1) n r(n+oc+1)/n! r(a+1) , 

where r denotes the Euler's gamma-function, the Legendre polynomial L n coincides with 
J° , while the Chebyshev polynomial T n is equal to 4 n [(n!) 2 / (2n)!] J~ 1/2 . Finally, we 
recall that the J“ , n > 0, are the eigenfunctions of a Sturm-Liouville operator, more 
precisely they satisfy 

(11.1) ((1-C 2 )“ +1 J*')' + n(n+2a + 1) (l-C 2 )“J“nO . 

We refer to [DR, § 1 . 1 3] for the properties of these orthogonal polynomials. 

Next, let N be a fixed integer > 3. We denote by C* , 0 < j < N, the zeros of the 
polynomial ( 1 -C 2 ) Jjj'. with -1 = Co < £?<•■•< = 1 • There exist weights , 

0 sj j < N, such that the Gauss- Lobatto quadrature formula 

(11.2) /.‘^(C) ( 1 -C 2 )“ dt; - Z^j N =0 ) ef 
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is exact for any polynomial in P 2 n_i(A). We shall need the interpolation operator vjj 
associated with these nodes : for any function q> in ($°(A), i^<p belongs to P N (A) and 
satisfies 

(11.3) v A <pU A ) = «p(t A ) ,0<j<N . 


Remark 11.1 : It is well-known that the zeros cj" and the weights pj" satisfy 

(11.4) [ Cj C = cos((N-j)n/N) , 0 < j < N , 

pj" = n/N , 1 < j < N- 1 , and p£ = p£ = n/2N 
Note that, although the Arccos of the tj . 0 < j < N, are not strictly equidistant, their 
distribution is coarsely the same. 

Example 11.1 : For N = 15, we give the values of - t A = , 1 < j < 7. 



Legendre case 

Chebyshev case 

j- 1 

0.9695680462702179 

0.9781476007338056 

) = 2 

0.8992005330934721 

0.9135454576426009 

j = 3 

0.7920082918618152 

0.8090169943749474 

j -4 

0.6523887028824931 

0.6691306063588582 

j = 5 

0.4860594218871376 

0.5000000000000000 

j = 5 

0.2998304689007632 

0.3090169943749474 

J = 7 

0.1013262735219494 

0.1045284632676535 


Now, let us consider the two-dimensional domain 0 = )-1 ,1 [ 2 . For any nonnegative 
integer n, we denote by P n (G) the space of restrictions to Q of polynomials of degree < n 
with respect to each variable, i.e. the space P n (A) <s> P n (A); we also introduce the space 
P*(0) of polynomials of P n (Q) which are equal to 0 on the boundary 90. 


For the fixed integer N, we define the grid ZjJ by 

(11.5) Z A = { x jk A = (cf ,t k A ) ;0< j,k<N} . 

The idea of defining the grid from the nodes of a Gauss type quadrature formula was first 
presented in [Go]. 

To each point x jk in , we associate the weight pj A = p A p A . That allows us to define 
the following bilinear form on 6°(0) x C°(0) 

(11.6) («p,i|/) a . n = Lj N i0 £k=o Pjk • 

Since the quadrature formula (11.2) is exact on P 2 N _i(A), it coincides with the scalar 
product (,,.) A on P N _,(0); it is known [CQI , §3] that, on P N (0), it is still a scalar 
product, and the norm : ip -► (tp.ip)^ is equivalent on P N (0) to the norm : <p ->■ ||(p|| 0 Afl , 
with equivalence constants independent of N. Finally, we define the interpolation operator 
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3^ on the grid z£ in the following wag : for any function <p in 6°(0), 3,Jip belongs to 

P N (0) and satisfies 

(11.7) 3^<p(x) = tp(x) ,x€ Z„ . 


11.2. Th e c ollocation prQ b1 m .f.QrJiomt>geneQus boundar y c on dition?. 


We choose the space X N of discrete velocities equal to P N (0) 2 and the space M N of 
discrete pressures equal to a subspace of P N (0) which we will precise later. 

Let us assume that the data I belong to C°(Q) 2 . The collocation approximation of the 
homogeneous Stokes problem (1. 3)0.4) is the following : Find a pair (u N , p N ) in X N x M N 
satisfying 


x c ~ A 
> * ^ “N 


(11.8) A (- v Au n + grad p N )(x) = f(x) , x € Z(J (1 0 

(div u N )(x) = 0 

together with the boundary conditions 

(11.9) A u n (x) = 0 , x € z£ 0 00 


In order to state a variational formulation of problem (ll.8) A (ll.9) A , we define 
three bilinear forms respectively on 6 2 (0) 2 x 6°(0) 2 , on C°(0) 2 x C'(0) and on 
O’(0) 2 x e°(0) by 

(11.10) V u € 6 2 (0) 2 , V v € C°(0) 2 , a AN (u,v) = - v (Au,v) an , 

(11.11) V v € 6°(0) 2 , V q € G'(0), b 1A N (v,q) = (v , grad q) A N , 

(11.12) V v <= ©'(O) 2 , V q<= C°(0), b 2AN (v,q) = - (div v , q) AN . 

Since the quadrature formula (11.2) is exact on polynomials of degree < 2N- 1 , we have 
(11.1 3 ) l V u € P N (0) 2 , V v € P^(0) 2 , a LN (u,v) = v (grad u , grad v) LN , 

( 1 1 . 1 3) c V u € P N (0) 2 , V V € p„(0) 2 , a C N (u,v) = v (grad u , grad (vw) w _1 ) A N 

Moreover, we note that, in the Legendre case, the two forms b 1LN and b 2LN coincide on 
P^(0) 2 x P N (0) while, in theChebyshev case, one has 
(11.14) V v € P|J(G) 2 , V q € P N (0), b 1CN (v,q) = - (div (vw) w _1 ,q) CN . 


Now, it is clear that problem (ll.8) A (ll.9) A is equivalent to the following variational 

A 

one: Find a pair (u N , p N ) in P^(0) x M N such that 


(11.15). 


V v N € P^(0) 2 , a AN (u N , v N ) + b 1A N (v N , p N ) = (f,v N ) A N 
V<*€P M (0), * 3 2A,n^ u n • 9m) = 9 


Finally, to discretize the Navier-Stokes equations (1. 0(1.4), let us consider the 
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nonlinear term. Since the solution u of (1. 1 ) is divergence-free, the quantity (u. V )u is 
equal to £j =1 3(u i u)/5x i , where u, and u 2 are the components of the velocity u. Though this 
property is no longer true for the discrete problem in the general case, it seems more 
convenient to choose the second form, for reasons of numerical stability. Moreover, if a 
continuous function u is known by its values at the nodes x of Z^ , it is easy to derive the 
values of u ( u at the same nodes, hence to compute 3 l J‘(u i u). The pseudo-spectral 
approximation, as suggested in [01], consists in differentiating this interpolant i.e. , in 
replacing dCUjiO/dXj by 93£(UjU)/9Xj . 

These two arguments lead us to the following discretization of the Navier-Stokes 
equations : Find a pair (u N , p N ) in X N x M N satisfying 

(11 1 6) a [- v Au n + grad p N + L, 2 =1 93^(u Nj u N )/9Xj](x) = f(x) , x € Z(J Pi 0 , 

(div u N )(x) = 0 , x€Z„ , 

together with the boundary conditions (ll.9) A . 

Of course, this system is equivalent to the variational problem : Find a pair (u N , p N ) 
in P^(Q ) 2 x M n such that 

V v N € P,;(0) 2 , a AN (u N , v N ) + b 1AN (v N , p N ) 

(II.17) a I + H-j (03 N (u N j u N )/9Xj , v n ) a _n = (** V N^A,N ’ 

V q N € P N (0), b 2A N (u N » q^,) = 


Our purpose is to choose appropriate discrete pressure spaces M N , such that 
problems (ll.8) A (ll.9) A and (II. 1 6) A (ll.9) A are well-posed. 


11.3. The discrete pressure spaces. 

It is known [Mo][Me][BMM][BCM] that the space P N (0) contains "spurious" modes 
for the pressure, i.e. , polynomials the gradient of which vanishes at the collocation nodes of 
Z^ f| 0; of course, even if they can be solved numerically (see Section V), the collocation 
problems cannot be well-posed if any of these modes belongs to M N , hence we have to 
characterize them. More precisely, for i = I or 2, we define the subspaces 

(11.18) Z iA N = { q N € P N (0); V v N € P^(fi) . b iA N ( y N , q N ) = 0 } 

Let us also introduce, in the Chebyshev case, the polynomial s N of P N ( A) which satisfies 

(11.19) V«p N €P N (A), s N (tf ) <P N (tj C ) $[ = /', <P N (C)dt . 

Finally, we need the Lagrange polynomial r A associated with each node Cj A , 0 < j < N : rj A 
belongs to P N (A), is equal to 1 in t A and to 0 in any other node C A .0<ksgN,k*j. 
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We have [BMM , Lemma V. 1 ][BCM, Prop. V.2 and Y.3] 

Lemma 11.1 : The sub space Z iAN , ini or 2 , is of dimension 8. It is spanned 

1 ) in the Legendre case, for i = 1 or 2 , by { L 0 , L N }® 2 and { rj; , r{j }® 2 ; 

2) in the Chebyshev case, for i = 1 , by { T 0 , T N }® 2 and { r£ , rj$ }® 2 , for i = 2 , by 
{ s N , T n }® 2 and { r c 0 , r£ }® 2 . 

Let M iA N denote the orthogonal subspace of Z iAN in P N (0) with respect to the scalar 
product (,,.) AN • In the sequel, we shall always choose the space of discrete pressures M N 
such that the orthogonal projection operator n N : M N -+ M 1AN with respect to the scalar 
product (.,.) AN satisfies 

(11.20) V q^j € M N , Ildfjllo.A.n ^ c H n n ^ o.A.n • 

Remark 11.2 : Of course, the choice M N = h 1A N is the most natural one. However, this 
space has not good approximation properties since it can be checked that all its elements 
vanish in the corners of the domain 0 (which is a priori not the case for the exact 
pressure). On the opposite, for any real number X, 0 < X < 1 , it is possible to build 
subspaces M N which satisfy (11.20) and such that the following inclusion holds 

(11.21) P [xn] (0)cM n ( [XN] denotes the integral part of XN) , 

which implies that these M N have good approximation properties. Examples of such spaces 
are given in [BMM , Prop. V.3] in the Legendre case and in [BCM , (IV. 6 1 ) and (IV. 49)] in 
the Chebyshev case. 


Now, problem (ll.8) A (ll.9) A seems overspecified, since there are eight equations 
more than unknowns. Due to the definition (11.18) of Z 2A N * it turns out that. problem 
( 1 1 . 1 5 ) A is equivalent to : Find a pair (u N , p N ) in P^(0 ) 2 x M H such that 

(11.22) A | Vv n €P*( 0) 2 , a A . N (u N ,v N ) + b 1AN (v N ,p N ) = (f,v N ) AN , 

^ Qn € ^2A,N ■ ^2A.N^ U N • = ® 

Clearly, the continuity equation in (ll.8) A is redundant. However, let us denote by S e the set 
of the four corners of the square 0, and introduce a set S A of four collocation points in 
\ $ c satisfying the following property : 

(11.23) det ( q K (Xj) ) * 0 , I<J,K<S4 , 

where x 0 runs through $ A and q K runs through { L 0 , L N }® 2 in the Legendre case and through 
{ Sn , t n }® 2 in the Chebyshev case. The following result is proven in [BMM , Prop. V. 1 ] and 
[BCM, Prop. V.7], 
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Proposition 11.1 -.Assume that hypothesis (11.23) holds. Problem (ll.8) A (M-9) A is 
equivalent to the following one : Find (u N , p N ) in X N x M N satisfying 
( 1 1 . 2-4) A (- v Au N + grad p N )(x) = f(x) , x € z£ D 0 , 

(div u N )(x) = 0 , x € ZjJ \ { S c U S A } . 

together with the boundary conditions (ll.9) A . 


In the same way, we have 

Proposition 11.2 : Assume that hypothesis (11.23) holds. Problem ( 1 1 . 1 6) A ( 1 1 . 9) A is 
equivalent to the following one : Find (u N , p N ) in X N x M N satisfying 

(11.25). 


[- v Au n + grad p N + £. :1 a3^(u Ni u N )/ax ( ](x) = f(x) 
(div u N )(x) = 0 , x € Z^ \ { $ c U S A } . 

together with the boundary conditions ( 1 1 . 9) A . 


In Section V, the reader will find practical ideas for solving this system, in particular 
how to choose a convenient set of degrees of freedom for the pressure. 


11.4. The collocation problem for inhomogeneous boundaru conditions . 

In this paragraph we assume that f belongs to 6°(0) 2 and that the boundary data g are 
such that the gj = g| fj , J € Z/4Z, satisfy (1. 1 5) and (1. 1 6) but also belong to C°(Fj) 2 and 
satisfy 

(11.26) gj(aj) = g j+1 (aj) ,J€Z/4Z . 

We are now interested in the approximation of problem (1.3)0. 2). The first idea is to use 
the same discrete problem (ll.8) A as in the homogeneous case and simply replace the 
boundary equation (ll.9) A by 

(11.27 ) A u N (x) = gj (x) ,X€ Z A nr, , J € TL/ATL . 

But it turns out that this problem has no solution in the general case. Indeed, if the equation 
div u N was satisfied in any point of Z^ , we would derive div u N = 0 exactly. In particular 
this would imply five conditions for u N at the boundary : 

(11.28) (div u N )(aj) = 0 ,J€Z/4Z , 

(11.29) Z.j € z/ 42 jfj U N ■ n j da = 0 . 

These equations solely depend upon the values of u N at the boundary, hence upon the 6jJgj , 
J € TL/ATL. In general they are not verified, even if (1. 1 6) holds (examples of functions gj 
satisfying our assumptions but violating (11.28) and (11.29) are given in [BCM, (Y.6)] and 
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in [Me, Chap. 4]). 

That is why we propose the following discrete problem : Find (u N , p N ) in X N x M N 
satisfying (ll.22) A together with the boundary conditions (ll.27) A . 

Note however that this last problem is not so far from a collocation one, as the 
following proposition states it. 

Proposition 11.3 : Any solution (u N , p N ) of problem (ll.22) A (ll.27) A in X N x M N 
satisfies the collocation equation 

( 1 1 . 30) A (-v Au N + grad p N )(x) = f(x) ,x(Z N A f)0 . 

Remark 11.3 : By noting that the space M 2A N is exactly the image of P^(Q) 2 by the 
divergence operator, it can be seen that solving the equation 
V q € M 2A N , b 2A N (u N , q) = 0 

in (1 1 . 2 2 ) A is actually equivalent to the minimization of ||div u N || AN ; this condition is 
implemented in practice, as will be seen in Section V (cf. also [Me]). 


(11.31), 


In the same way, we define the approximation of the inhomogeneous Navier-Stokes 
equations (1. 1 )( 1 .2) as : Find (u N , p N ) in X N x M N such that 

V V N € P|J(0) 2 | a A,N^ U N * V N^ + ^1A,N^ V N ■ Pn^ 

+ 2--i =1 (93 N (u Nj U N )/dX. , V N ) A N = (f,V N ) A N , 

V € ^2A,N > b 2A.N^ U N > 'W = 0 

together with the boundary conditions (ll.27) A . We also have the 

Proposition 11.4 : Any solution (u N , p H ) of problem (11.3 1 ) A (ll.27) A in X N x M N 

satisfies the collocation equation 

(ll.32) A [-v Au n + grad p N + Lf =1 (93^(u Ni u N )/aXj](x) = f(x) ,xCZ N A flfi 


In the following sections, it will be proven that the four discrete problems are 


well -posed. 
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III. Convergence results for the Stokes problem. 


The convergence of the method in the case of homogeneous boundary conditions has 
already been thoroughly analysed [BMM][BCM], hence we only recall the results. Then, we 
extend them to the nonhomogeneous case. 


IM.l. The case of homogeneous boundaru conditions. 

Problem ( 1 1 . 8 ) A (ll.9) A is actually analysed through its variational formulation 
( 1 1 .22 ) a . We recall the main properties of the bilinear forms involved in this formulation, 
which are the corner-stone of the study. For i = 1 or 2, let us define the kernels 
(HI O <iA.N = { v n ^ Pn((0 2 ; V q N € P N (0), b iAN (v N , q N ) = 0 } . 

Clearly, K, A N in the Legendre case and K 2A N in both cases coincide with the subspace of 
divergence-free polynomials in P^(0) 2 , while K, A N in the Chebyshev case is the subspace 
of polynomials v N in P^(Q) 2 such that div (v N u) is equal to 0. 


(III. 2) 


Proposition III. 1 : There exist constants * , 6, and 6 2 independent of N such that 
forms a AN , b 1AN and b 2A N satisfy the following continuity properties 

V U N e P N (0) , V V N €P N (0) , I a A ,N ( U N > V N ^ I ^ If II U N I1 1 ,A ,0 H V nHi,A,0 

^ v n € P*n(^0 , v q N €P N (Q), 0 1A N (v N , q N )| ^ 8, ||v N II, A n llq N ll 0 ,A.n 

^ v n € ^ ‘-In € P N ( 0 ) , |b 2A N (v N , q N )| < 8 2 ||v N || 1A n ||q N || 0An 

In the Legendre case, there exists a constant a L > 0 independent of N such that 
( 1 1 1 • 3)|_ ^ W N € P N (0) , a UN (w N , w^) > oc L || || ^ ^ ; 

in the Chebyshev case, there exists a constant oc c > 0 independent of N such that 

(1 1 1 . 3 ) c v W N € K 2A n , 3 V N 6 K 1A n . V N ^ 0 / 9 C,n( W N > V N^ ^ a c II W N lll,u,0 II V N Ht.u.O 
There exists a constant §,> 0 independent of N such that 

(111. 4) Vq N € M^ N ,3v N € P’(Q) 2 ,v N *0/ 

biA.N^N ’ C W ^ N H V N II 1 ,A.O H^N Ho.A.O 


the 


if hypothesis (11.20) holds, there exists a constant ^ > 0 independent of N such that 

(111.5) V q N € M N , 3v N € P‘(Q) 2 , v N * 0/ 

Vn( v n ’ ^ M II v n I1 1 ,a.o II c In Ho.A.n 


Using a well-known theorem for saddle-point problems [B ] [ GR , Chapter I , Corollary 
4.1 ][T, Chap. I, Th. 2.1] in the Legendre case and its generalization to nonsymmetric 
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problems [N][BCM, Corollary 11.2], we derive from this proposition the following 
theorems [BMM, Thm V. 1 ][BCM,Thm V. 1 ]. 

Theorem III. 1 : Assume that hypothesis (11.20) holds. For any function f in 6°(0) 2 , 
the collocation approximation (ll.8) A (ll.9) A to the Stokes problem (1.3)0. 4) has a 
unique solution (u N , p N ) in P^(O) 2 x M N . 

The error estimates have been proven respectively in [BMM, Thm V. 1 and V.2] in the 
Legendre case and in [BCM, Thm V.2 and V.3] in the Chebyshev case. Note that the main 
arguments of the proofs will be recalled in the following subsection, in order to study the 
inhomogeneous case. 

Theorem 1 1 1.2 : Assume that hypothesis (11.20) holds, that the solution (u,p) of the 
Stokes problem (1. 3)0.4) is such that u belongs to H A (0) 2 for a real number s^ 1 , 
and the data f belong to H A (Q) 2 for a real number o > 1 . Then the solution (u N , p N ) of 
problem (ll.8) A (ll.9) A satisfies 

(1 1 1.6) llu-uj, A 0 < c ( N 1 -* Hull, A0 + N 1 * 2 *-' ||f || ffAft ) 
for a constant c independent of N. 


Theorem II 1.3 : Assume that hypotheses (11.20) and (11.21) hold and that the solution 
(u,p) of the Stokes problem (1. 3)0.4) belongs to H A (Q) 2 x H^'^O) for a real number 
$> 1 , and the data f belong to H A (G) 2 for a real number a > 1. Then the solution 
(u N , p N ) of problem (ll.8) A (ll.9) A satisfies 

(III. 7) IIP-PnIIo,A, 0 ^ C { N 3 * ( ll u ll$,A,0 + IIpII*- 1,A,0 ) + N * " ll*ll<r,A,0 ) 

for a constant c independent of N. 


Remark III. 1 : Let us consider for a while the problem : Find (u N , p N ) in P^(0) 2 x M N 
such that 

( 1 1 1 . 8 ) A V V N € P^(o) 2 , 3 a.N^ U N > V N^ + ^1A,N^ V N 1 Pn^ = ^*> V n)a > 

v q N € P N (0), t>2A,N^ U N > = ® 1 

which is problem (II. 15) A with 0,v N ) AN replaced by 0,v N ) A . Then, Theorems III. 1 to III. 3 
are still valid. Furthermore, by reading the proof of [BMM][BCM], it is easy to see that the 
estimates (II 1.6) and (II 1.7) can be replaced respectively by 


011.9) 

and 


I u_u nIIi,a,o ^ c N 


i-»i 


lull 


s.A,0 


011.10) Hp-pJIoao < c N 3 ‘ s ( llu||,,A.o + IIPll,-1.A.ft ) 
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This will be used in the following section. 

1 1 1.2. The case of inhomogeneous boundaru conditions. 

Our purpose is now to study the discrete problem (ll.22) A (ll.27) A . Since we need an 
element in the space of trial functions X N that satisfies, in a discrete sense, the boundary 
condition (1.2), we state the following lifting result that can be derived from [BM1 , 
Prop. V. 1 ]. 

Lemma 111.1 : There exists an operator QjJ from the subspace of all polynomials 
= ^ ^nj )j«z/ 4 z ’ n Hjszmz ^n^j) satisfying 

(111.11) tpfjj(3j) = <pf| j + j (3j) , J € 2/42 , 

into P N (0) such that, for any such polynomial $ N , 

(111.12) Q^n) = <P NJ on Tj , J € 2Z/4Z . 

Moreover the following estimate is satisfied 

(111.13) II<WII 1 .a.o < c N ' “^j«z/4zll , * , Njllo.A.r J 

Sketch of proof : From [BM1 , Prop. V.l], we derive that there exists an operator QjJ 
which satisfies (III. 1 2) and such that one has 

HQn^n^Iiao ^ c (N Z^jez/4z HfNjllo,A.rj + ^ 2 " ^j€z/4z ^ 

Then, estimate (111.13) follows from the previous line and the inverse inequality [Q, (2.4) 
and (3.2)], valid for any polynomial <p N in P N (A), 

H'PnIIl^a) <cN ’ “H'PnIIo.a.a • 

The previous result allows us to check that the discrete problem (ll.22) A (ll.27) A is 
well- posed. 

Theorem 1 1 1.4 : Assume that hypothesis (11.20) holds. For any function f in C°(Q) 2 , 
the collocation approximation (II. 22) A (11.27 ) A of the Stokes problem (1.3)0. 2) has a 
unique solution (u N , p N ) in X N x M N . 

Proof : If we choose u N b equal to the image of ( i^gj )j €2/42 b y the operator , the 
polynomial u N = u N - u Nb belongs to P^(G) 2 . Then the pair (u N , p N ) is a solution of 
(ll.22) A (II.27 ) a if and only if the pair (u N , p N ) satisfies : 
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V v N € P^(0) 2 , a AN (u N , v N ) + b 1A N (v N , p N ) 

<"'.14)a = ( <>V nVn ~ a A,N^ U N,b • V N^ > 

V € M 2A,N > b 2A,N^N - ^ = " b 2A,N^ U N,b ’ ' 

Due to Proposition III. 1 , we derive the result from [BCM, Corollary 1 1.2]. 

Next, we will study the approximation of divergence-free functions by divergence- 
free polynomials, thus generalizing the result of [SV] to the case of non homogeneous 
boundary conditions. Let us set 

(111.15) K a (G) = { w € H A (0 ) 2 ; div w = 0 in 0 } . 

Lemma III. 2 : There exists an operator from K A (Q) into P N (0 ) 2 D K A (Q) such 
that the following estimate is satisfied for any real number s > 3 : for any function w in 
K a ( 0 ) D H a ( 0) 2 , 

(111.16) llw-Rj^wll, A 0 < c N 1 "*||w || sA0 . 

Proof : Let us recall [M I , Remark 11.3 and Lemma IV. 2] that there exists an operator tt ^' 2 
from H a (A) into P N (A) which satisfies for any function ip in H a (A) 

(111.17) n(J' 2 ip(- 1 ) = <p(-1 ) and n£ , 2 ip( + 1 ) = <p( + 1 ) , 

(dn^' 2 4>/dC)(- 1 ) = (d<p/dt)(-l) and (dn^- 2 ip/dC)( + 1 ) = (d«p/dC)( + 1 ) . 
Moreover, we have for any function <p in H a (A), s > 2, 

( 111 . 18 ) |j «P -nj^' 2 tp H 2 A A + N il‘P -n N ' 2 t Plll,A,A + N 2 H l P -n N' 2 | Pllo,A,A C N 2 ' S ||<p || $ A A . 
Next, for any function w in K A (Q), there exists a unique function <p in 

H a (Q) fl L ao ( 0) such that w is equal to curl 41 in Q; moreover, if w belongs to H A ( 0 ) 2 , 
s s* 3, it satisfies 

(111.19) lkll st 1 A 0 <c||w || sA0 . 

Setting tp N = (n A, 2 ®njJ , 2 )«p, we define R,Jw as being equal to curl <p N . It remains to estimate 
II W— R n w||, A ,fj < c ll'p - ^ 2 ,A.O 

<c( ||ip-ip N || H 2 (AL 2 (A)) + lli , -'i' N ll H |( AiH ^( A )) + ll4'-q' N ll L |<A ,h a (a)> )• 
From (111.18), we infer 
lk-(n N A ' 2 ®n N A ' 2 ) 4 <|l H 2 (AL 2 (A)) 

< ll«P-(TT^- 2 <S>td)v|/ || h 2(a ,l|(A)) + ll n N 2 ®^ id-rl N ^H a (A ,L a (A)) 
^cN 1 ‘ s |kl| H s + i (AL 2 ( A ) ) + ||id<E>(id-n^' 2 )i|»|| H 2 ^ A>L 2 ^ A ^ 

<CN 1_S ( ||ip|| H s + 1 (AiL 2 (A)) + ||>p|| H 2 (A H s-1 (A)) ) 
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< C N 1-3 ( ll*l H ^ 1 ( AiL *(A)) + M H a (A,H a _1 (A)) > 

<CN ll^Ust I.A.ft • 

The quantities lk-*J H ’(A,Hj(A» and l! L 2( A .h|(A)) are estimated in the same way. 
Finally, we obtain 

J| w— R n w ii i , a,o ll*i* II»+ i,a,o ■ 

which together with (111.19) gives the result. 

Of course, other divergence-free polynomials approximations of divergence- free 
functions in H A (0) 2 can be built, However, note that, for any function w in K A (Q), the 
operator satisfies 

(111.20) (RjJwXaj) = w(aj) ,J€2/4Z . 

Moreover, it has the following useful property. 

enrol lari , I 111.1 : The operator RjJ satisfies the following estimate for any real number 
x ^ 2 : for any function w in K A (G) such that the trace W| fj belongs to H A (fj) 2 , 

(111.21) ll w (rj“ R N w lrjHo.A.rj ^ c N T H w |rjllt.A.rj 

Proof : We write 

||w| rj -R£ w |r J 0 Afj < c ( ||dq;/ax-aq/ N /dx|| 0 A rj + ||0ii//ay-a<|/ N /ay|| OArj ) . 

In the case J = I or III for instance, using (ill. 1 7) and noting that the operators d/dx and 
id®njj' 2 commute, we have 

(a(n^' 2 ®njj ,2 )^/9x)(± 1 ,y) = (a(id®n£ -2 )«|i/ax)(+ 1 ,y) = (id®n|J' 2 )(a<ti/0x)(± 1 ,y) 
and, similarly, 

(a(n£ ,2 ®n^- 2 )i|</dy)(± 1 ,y) = (a(id®n^ ,2 )i|i/ay)(± 1 ,y) , 

so that, by (111.18), 

H w |rj -R N w |rj Ho,ajj < c < «(3t/ax)-(id®n N A - 2 )(a4./ax)|| 0A _ rj 

+ ||<|i-(id®njj- 2 )i|i|| ) Arj ) 

< c n* t ( ||a^/ax|| T A rj + M t+1Arj ) < c n- t llw| rj || tArj . 

The cases J = II and IV are studied in the same way by exchanging the variables x and y. 

We are going to introduce a slightly different approximation to the Stokes problem 
(1. 3)0.2) , that satisfies conditions (1.16) and (11.26). In all that follows, we assume that, 
if (u,p) is the solution of (1. 3)0.2), the function u belongs to H A (Q) 2 ; then, we set 

(111.22) z N . t = R N A .,u . 
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Let us remark that the data z. 


N-l|fj 


, J € TL/ATL, define a boundary condition that satisfies 


(1.16) and (11.26). Hence we can define two auxiliary problems 
a) Find (u,p) in H A (0) 2 x L A0 (Q) such that 
( 1 1 1 .23) A Vv(H| 0 (fl) 2 , a A (u,v) + (v , grad p) A n (f,v) A , 
V q € L a (Q), (div u , q) A = 0 , 

and satisfying the boundary conditions 


(111.24)* u = z, 


N-i|rj 


on r | , J € Z/4ZL 


b) Find (u N , p N ) in P N (0) 2 x M N such that 

(lll.25) A V V N € P|^(0) 2 , 3a,N^N • V N^ + ^1A,N^ V N 1 Pn^ = ^’ V N^A,N > 

^ 9n € ^2A,N ’ b 2A.N^N • %) = 0 . 

and satisfying the boundary conditions 

(III 26) a u N (x) = z N _,| rj (x) , x € f) l"j , J € TL/ATL . 

The error bound between the solutions u and u N of problem (1.3)0. 2) and (ll.22) A (ll.27) A 
will be obtained by studying the differences between u and u, D and u N , u N and u N . 


Lemma III. 3 : Assume that the solution (u,p) of the Stokes problem (1.3)0. 2) is such 
that u belongs to H A (Q) 2 for a real number s ^ 3. The following estimate is satisfied 
011.27) l|u-u||, A0 < c N 1 ' s ||u|| sA0 

for a constant c > 0 independent of N. 


Proof : Since the pair u-u is the solution of a Stokes problem with null body forces and 
boundary data equal to U| rj - * N _i| rj . J € TL/ATL , it follows from the stability estimate 
0.22) that 

H u_ “lli,A,n ^ c £j€Z/4z ll u |rj -z N-i|rjll(i-«)/2,A,rj • 

Due to the trace theorem [LM, Chap. 1 , Th. 8.3][BM 1 , Thm 11.2], that implies 

H u - D IIiao< c II u - z n-iIIiao • 

Then, we deduce the lemma from Lemma 1 1 1.2. 

Similarly we can obtain an error bound between u N and u N . 

Lemma III. 4 : Assume that the boundary data gj , J € TL/ATL, belong to H*(l" j) 2 for a 
real number x > 2. The following estimate is satisfied 
(111.28) II u n _ u n |Ii A ,n < c N 7/2 T ZI-jez/42 HAjIItAIj 
for a constant c > 0 independent of N. 
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Proof : It follows from problems ( 1 1 . 2 2 ) A (1 1 . 2 7 ) A and ( 1 1 1 . 2 5 ) A ( 1 1 1 . 26 ) A that the 
polynomial u N -u N is the collocation approximation of a Stokes problem with null body 
forces and boundary data equal to ^ 6 2Z/4Z. Let w N b denote the image of 

< 9j - z N-i|rj )j e z/ 4 zr b y the operator QjJ . Setting w N = u N - 0 N - w MJb and r N = p N - p N , 
we see that the pair (w N , r N ) is the only solution in P^(0) 2 x M N of 
(IN 29) a | VV N (P*(Q) 2 , 0 A .n( W N -V + W V N ■ = - a A.N (w N.b • V N ) ■ 

V q N € M 2A n , b 2A N (w N , q N ) = - b 2 A ,j/ w N.b * 

Using [BCM , Corollary 11.2] together with Proposition 1 1 1. 1 , we obtain 

II W N I1 1 , A ,0 ^ C N II W N,b 111 , A ,0 > 

so that, using Lemma 1 1 1.1 , 

II u n“^n lli, A ,n ^ c N 2 N 1 ( llflj— Sj Ho.A.rj + H u |rj~ z N-i|rj Ho. A ,rj ) • 

The lemma follows from Corollary 111.1 and from the following estimate for the 
interpolation error [CQ I Thms 3. 1 and 3.2] , valid for any real number s > 1 /2 : 

(in.30) II<p-<*pII 0A a <cN ,/2+ “- s ||<pII sAA • 


Finally, in order to get now an error bound between u and u N , we note that problem 
( 1 1 1 . 25) A ( 1 1 1 . 2 6) A is a discrete approximation of problem ( 1 1 1 . 2 3) A ( 1 1 1 . 2 -4) A . That 
allows us to derive the following estimate. 

Lemma 1 1 1.5 : Assume that the solution (u,p) of the Stokes problem ( 1 . 3)( 1 .2) belongs 
to H a ( 0) 2 x H a _, (0) for a real number s ^ 3 , that the data f belong to H A (Q) 2 for a 
real number o > 1 and that the boundary data gj , J € Z/4Z, belong to H A ( r j ) 2 for a 
real number x ^ 2. The following estimate is satisfied 

(111.31) l|u-u N II, Afi < c ( N 1 5 l|u|| s A 0 + N 2 “ ||f|| ffA0 + N T L Je2/42 llgjll tArj ) 

for a constant c > 0 independent of N. 


Proof : Let us set u* = u - z N _, and = u N - z N _, . Thus, (u*,p) is the solution in 
H; 0 (O) 2 X L 2 0 (O) of 

( 1 1 1 . 32) A | VveH>) 2 , a A (u*,v) + (v , grad p) A = (f,v) A - a A (z N _, , v) , 

V q € L A (Q), (div u* , q) A = 0 , 

and (ujj , p N ) is the solution in P^(G) 2 x M N of 
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( 1 1 1.33), 


^ V N € Pn(«) ’ a A.N^ U N ' V N^ + ^1A,N^ V N > Pn^ ~ ^> V N^A,N " a A.>/ Z N-1 ’ V N^ • 


1 V q N € M 2A N , b 2A N (ujJ , q N ) _ 0 . 

Next, we use an abstract error estimate due to [BCM, Corollary 11.3] : since the forms 
a A (.,.) and a A N (.,.) coincide on P N _,(0) 2 x P *(Q) 2 and since z N _, belongs to P N _,(0) 2 , for 
any polynomial v N _, in P^,(0) 2 and any polynomial w N in P„(Q) 2 such that 
V q N € M 2A n , b 2A N (w N , q N ) = 0 , 

the following estimate holds (note that u* - uJJ = u - u N ) 


(111.34) 


l“-Ml.A.O< C <l|u‘-V N _ 1 || )Aft + ||u*-W N || 1A0 

(* > " (* ’ z n\,n 


sup 

z n € Pm(O) 


N' 


C N AO 


The more convenient here is to choose v N _, = w N = 0. So it remains to estimate ||u* ||, A 0 and 
the last term. 

I ) We have 


ll U *lll,A.O< l|u-u||,,A.O + ll U ' Z N-llllAO • 
so that, by Lemmas 1 1 1.2 and 1 1 1.3, 

(111.35) llu‘11, Aft < C ( N 1 4 l|u|| sA0 + N T 2 Lj € 2/42 HSj Ht.AXj ^ 

2) We recall [CQ1 , §3] that the scalar product (.,.) A induces a norm on P N (0) which is 
equivalent to ||.|| 0 A0 . Hence, choosing f N _, in P N _,(0) 2 , we obtain for any z N in P„(0) 2 

- Z N^A “ ’ Z N^A,N = - *N-I > ^A “ * “ *N- 1 ’ Z N^A.N 

< C ( ||f - f N _ -i Ho AO + II* ~ * Ho,A.O ^ Ho AO 

Let us recall that the orthogonal projection njjfrom L 2 (Q) onto P N (0) satisfies the 

following estimate for any <p in H A (Q), s > 0, 

(Hi. 36) H<P-n^|| OAft ^cN- s ||«p|| sAO . 

Taking for instance f N _, = TT , f and noting that 3^ is equal to <s> , we derive from 


(111.30) and (111.36) 

(111.37) (» ,z N ) A - (I ,z N ) AN <cN u2ot ' <r 


H*H<r,A,0 H^Ho.A.O 


Finally, estimate (111.31 ) follows from (111.34), (111.35) and (111.37). 


From Lemmas III. 3 to III. 5, we derive the main error estimate. 

Theorem III. 5 : Assume that hypothesis (11.20) holds and that the solution (u,p) of the 
Stokes problem (1.3)0. 2) is such that u belongs to H A (Q) 2 for a real number s^3, 
that the data f belong to H A (Q) 2 for a real number o > 1 and that the boundary data gj , 
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J <E TL/ATL , belong to H^Tj) 2 for a real number x>2. Then, the solution (u N , p N ) of 

problem (ll.22) A (ll.27) A satisfies 

(111.38) ||u-u N || , A 0 < c ( N ’-* || u|| s A 0 + N ' * 2 *' ff ||f ||, A 0 

+ ^ 2 'J 6 Z/4Z I! Oj lU.A.rj ^ 

for a constant c > 0 independent of N. 


We conclude with an estimate for the pressure. 

Theorem 1 1 1.6 : Assume that hypotheses (11.20) and (11.21) hold and that the solution 
(u,p) of the Stokes problem (1.3)0. 2) belongs to H A (Q) 2 x H A ' 1 (0) for a real number 
s > 3, that the data f belong to H A (Q ) 2 for a real number o > 1 and that the boundary data 
gj , J € TL/ATL, belong to H a (Tj) 2 for a real number x>2. Then, the solution (u N , p N ) 
of problem (ll.22) A (ll.27) A satisfies 
(111.39) IIP-P N ll 0 .A.n N ( H u ll$,A,n + II pH*— i.a,o ^ 

for a constant c > 0 independent of N. 


* J € 2/ 42 llOjIliArj 


} 


Proof : Using the Inf-Sup condition of Lemma ill. 1 , we derive from (1.2 1 ) A and (ll.22) A 
that, for any q N in M N , 

Hp-Pn Ho.a. 0 ^ c ^ f II u u m Hi.a.o + Hp-^n Ho.a.o 

r > 9 rad ‘Wa ~ b 1A,N^ Z N ■ ^ f - Z N^A “ ^ ,>Z N^A.N ! \ 

+ SUp l jj jj + j; j] ] ) . 

Z N € P^(0) 2 H Z N Hl.A.O 

Owing to (11.21), taking for instance q N = n ( A XN] p and using (III. 36), (III. 38) and 
(111.37), we obtain easily (111.39). 


Remark 111,3 : As in the case of homogeneous boundary conditions, if we consider the 
problem : Find (u N , p N ) in X N x M N such that 

(lll.40) A V v n € P^(0) 2 , a AN (u N , v N ) + b 1A ^(v N , p N ) = (f,v N ) A , 

V q N € 0 2A N , b 2A N (u N , q N ) = 0 

and satisfying the boundary conditions (ll.27) A .Theorems 1 1 1 . 4 to I II . 6 are still valid, 
and it is easy to see that the estimates (111.38) and (111.39) can be replaced respectively by 


(111.41) Hu-uJI, Af) < c ( N s II u II s ,a,o + N £j€Z/ 4 Z HOj IU.A, rj ^ • 
and 

(111.42) l|P-P N llo.A.n { N ’ * ^ ll U lls.A.O + IIpIIs- 1^,0 ^ + N 11/2 T £.J € Z/4Z H 0J Ht.AXj ^ ' 
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IV. Convergence analusis for the Navier-Stokes equations. 


The aim of this section is to obtain, for the discrete problems (11.1 6) A ( II . 9 ) A and 
(li.31 ) A (ll.27) A , convergence results similar to those which were proven in the linear 
case. We begin by describing the main tools of the analysis, together with some properties 
of the exact equations. Then we establish some technical lemmas. This allows us to prove the 
convergence and to give error estimates for the velocity in both the homogeneous and 
inhomogeneous cases. Finally, error bounds are also derived for the pressure. 


IV. 1. The main tools 


To study the discrete problems, we shall use a fixed point theorem due to M. CROUZEIX 
[C, Th. 2.2], which is a refined form of the discrete implicit function theorem of [B RR ]. 
For the reader's convenience, let us recall this theorem : we consider a C-mapping F N 
from a Banach space Z N into itself and we assume that ujj is a point in Z N such that DF N (u^J) 
is an isomorphism of Z N . We denote by e N , x N and A n (t)), tq ^ 0, the quantities 


(IV. 1) 


e N - H F n( U nP Hz n *N - H(DF N (lljJ)) ll&(Z N .Z N ) 

A n (tj) = sup { ||DF n (w n )-DF n (u^) II^ (Zn . Zn ) ; W N € Z N and ||w N -u* || Z(J ^ t) } . 


Theorem IV. 1 : Let us assume that 2* N A N (2y N e N ) < 1 , then for each t ) > 2y N e N such 
that y N A n (tj) < 1, there exists a unique solution u N of the equation F N (u N ) = 0 in the ball 
$ N = { w N € Z N ; ||w N -ujJ || Zn t) }. This solution satisfies 

(IV. 2) V w N € S N , II u n -w n Hzn ^ 1 ~Xn A n (t)))] • H^n^ w n^ Hzn 


Let us precise in what framework we shall apply this theorem to the Navier-Stokes 
equations. We begin with the continuous problem. Let B A denote the subspace of all functions 
g in H A ‘“ >/2 (dO) 2 satisfying (1. 16) and (11.26). With the Stokes problem, we associate the 
operators T A and T A respectively from H A ’(0) 2 into H A0 (Q) 2 and from H A '(R) 2 x B A into 
H a ( 0) 2 : for any f in H A '(0) 2 , T A f is equal to the function u, where (u,p) is the solution of 
(1.3)0. 4) in H A0 (Q) 2 x L a0 (Q); for any (f,g) in H A '(0) 2 x B A , T A (f,g) is equal to the 
function u, where (u,p) is the solution of (1.3)0. 2) in H A (G) 2 x L A0 (Q). Of course, for 
any f in H A ’(Q) 2 , T A (f,0) coincides with T A f. 

Next, we consider the nonlinear term. We fix a function f in H^tQ) 2 and a function g 
in B a , and we define the following mappings 

(IV. 3) G(w) = a(WjW)/aXj - f and G(w) = ( G(w) , - g ) 


i 
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Clearly, the Navier-Stokes equations (I. 0(1. 4) have the following equivalent formulation : 
Find a function u in H A 0 (O ) 2 such that 
(IV.4) A u + T a G(u) = 0 . 

The Navier-Stokes equations (1.1 )( 1 .2) have the following equivalent formulation : Find a 
function u in H A (Q ) 2 such that 
( I V. 5 ) A u + T A G(u) = 0 . 

To check that these problems are well-posed, we need the 
Lemma IV. 1 : For any t in H A ’(Q) 2 , the mapping G is of class ($” from h’(0 ) 2 into 
H;’(0 ) 2 and from H^( 0) 2 into H A '(Q) 2 . Furthermore, for any w in H A (Q) 2 , the 
operator DG(w) is compact from H A (0 ) 2 into H A ’(0) 2 . 


Proof : Since the space H^fO) is contained in hVo), due to (1.7) and (1. 14), it suffices to 
prove that the mapping G is of class G°° from H’(O ) 2 into H^’(O ) 2 . For any u and w in 
h'( 0 ) 2 , we have 


V v € H^ 0 (O) 2 , H ,| 2 s1 } n (a(u i w)/ax.) vu dx| = | £. 2 =1 j 0 (u.w) 0(vu)/8x.) dx | . 
Since the mapping : <p -► ( | 0 |grad (ipu )| 2 u - ' dx ) ,/2 is a norm on H^fO) equivalent to 
the usual one [BM 1 , Lemma III. 2], we derive 

(IV. 6 ) V v<= H t J ) 0 (O) 2 , | Lf =1 /n (a(UiW)/ax.) vu dx | < c Lf =1 11^ w.|| 0ufl ||v|| 1uft . 
We recall [LM , Thm 4. 1 ][BM 1 , Lemma I II. 1 ] the imbedding of H 1 / 2 (0) into L 2 (Q), 
Moreover, using the Calderon extension theorem [A, Thm 4.32] together with [G, Thm 
1.4. 4. 2], we know that the mapping : (<p , 41 ) -»• cp is bilinear continuous from 
H 1 (Q ) x h'( 0) into H 1 _t (fl) for any e > 0. Hence we have for 0 < e < 1/2 
(IV. 7) llu. Wj || 0 w 0 < c ||u. Wj || 1/2 0 < c' luj Wj ||,_ e n < c" llu.H, 0 ||w j || 1 0 . 

From (IV. 6 ) and (IV. 7), we obtain 

V y €Hio(0) 2 . IL?., J 0 (a(u i w)/ax i )vudx|< c||u| 10 ||w|| t 0 ||v|| 1uft . 

Then, it is an easy matter to derive from (IV. 3) that G is of class C°° from H’(0 ) 2 into 
H‘'(0) 2 . The compactness of DG(w) from H 1 (0 ) 2 into H'^O) 2 , follows from the previous 
lines and from the compactness of the imbedding H ,_t ( 0 ) c H 1 / 2 (Q), 0 < e < 1 / 2 . 


Corollaru IV. 1 : For any 1 in H A '(0 ) 2 , problem ( 1 . 1 ) ( 1 . 4) has at least a solution 
(u,p )in H a 0 (O) 2 x L ao (0) . For any (f,g)/n H A '(0 ) 2 x B A , problem (1. 1)0.2) has at 
least a solution (u,p )in H A ( 0) 2 x L A 0 (O) . 
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Proof : In the Legendre case, the corollary states a well-known result [GR, Chapter IV, 
Thms 2.1 and 2.3]. Next, in the Chebyshev case, since the space H’’(G) is contained in 
H “ 1 ( 0 ) and the space B c is contained in B L , there exists at least a pair (u,p) in 
H 1 ( 0 ) 2 x Lg(G) solution of problem (1. 1 )(l.4) (resp. (1. 0(1.2)). From Lemma IV. 1 , G(u) 
is an element of H^’(G) 2 . Let (u',p') be the solution in H^(G) 2 x L 2 0 (G) of the Stokes 
problem with data - G(u). Then, both (u,p) and (u' ,p') are solutions in Hq(Q) 2 x L 2 (Q) of 
the Stokes problem with data - G(u); the uniqueness of the solution of the Stokes problem 
implies that u and u' coincide, and that p— p* is constant, equal to ( 1 /n) p(x) co(x) dx . We 
see that (u,p') belongs in fact to H^(G) 2 x L 2 0 (G) and is a solution of (1. 1 )(l.4) (resp. 
(M)(l. 2)). 

We state a last property of the continuous problem. It is interesting here to note that, 
since the second argument in G is constant, the operators 1 + T A DG(u) and 1 + T A DG(u) 
coincide on H A0 (G) 2 . 

Lemma IV.2 : For any real number q > 2, there exists a constant c(q,v) such that , if a 
solution (u,p) of problem (1. 1)0.4) (resp, (1. 1)0. 2 )) satisfies 
OV. 8) ||u|| Lq(0) < c(q,v) , 

the operator 1 +T A DG(u)/$ an isomorphism of H A0 (G) 2 (resp. the operator 
1 + T A DG(u))$an isomorphism of H A (G) 2 ). 

Proof : By the compactness result of Lemma IV. 1, the operator 1 + T A DG(u) is an 
isomorphism of H A0 (G) 2 and the operator 1 + T A DG(u) is an isomorphism of H A (G) 2 if 
and only if they are injective, i.e. the only solution (w,r) in H A0 (G) 2 x L A0 (G) of the 
following linearized Stokes problem 

V v 4 H A0 (G) 2 , a A (w,v) + (v , grad r) A + (DG(u).w , v) A = 0 , 

V q 4 L A (G) , (div w , q) A = 0 , 

is (0,0). In the Legendre case, the form a L is clearly elliptic on Hq(G) 2 ; in the Chebyshev 
case, it is proven [BCM , Prop. II 1.2] that, for any divergence-free function w in H^ 0 (G) 2 , 
there exists v in H^ 0 (G) 2 , satisfying div (vu) = 0, such that 
a^w.v) ^c||w|| 1un ||v|| l u0 . 

These properties, together with (IV. 6) in the Chebyshev case, give 
v H w ^ i , a,o ^ c =i II U i w j Ho.A.o 

Next, in the Legendre case, using the imbedding of H*(G) into any L S (G), s < +oo, we have 


at once 
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H u i w jllo.o< c (q> ll u ilUo ) ll w jlli,o ; 

in the Chebyshev case, a similar argument leads to 

II u i w jHo.u.n = II u > w j w ' /2 llo.n < c( q> llUjIUo) IIWj « 1/2 
In both cases, we obtain 


1 C (q) II Uj ll L q( 0 ) II W j I1 1 ,CJ ,0 


V ll W lll.A.O< c( q)ll U IUo)ll W lll.A.O • 

and the lemma is proven for an appropriate constant c(q,v). 


In the sequel, we shall always assume that the data f belong to a space H A (Q) 2 for a 
real number a > 1 and that the boundary data gj , J € TL/ATL , belong to a space H A (rj) 2 for 
a real number t > 2 and satisfy (1.16) and (11.26). We consider a solution u of the 
Navier-Stokes equations (I. I )( 1 . 4) (resp. ( 1 . 1 ) ( 1 .2) ) which is nonsingular in the following 
sense: the operator 1 + T A DG(u) is an isomorphism of H A0 (Q) 2 (resp. the operator 
1 + T a DG(u) is an isomorphism of H A (Q) 2 ); by virtue of Lemma IV. 2, such solutions exist 
for f and g small enough ! Even in the standard Sobolev spaces, regularity results of the 
solution (u,p) as a consequence of the regularity of f are not easy to derive [G, §7.3], 
whence we shall assume in the sequel that there exists a real number s ^ 1 (s > 3 in the 
case of non -homogeneous boundary conditions) such that the velocity u belongs to H A (0) 2 . 


We turn now to the discrete problems (II. 1 6) A (ll.9) A and (11.3 1 ) A (ll.27) A . As for 
the exact Navier-Stokes equations, we must define the operators T A N and t a N respectively 
from H a ’( 0) 2 into P^(O) 2 and from H A '(0) 2 x B A into P N (0) 2 : for any f in H A ’(0) 2 , 
T a N f is equal to the function u N , where (u N , p N ) is the solution of problem (II l.8) A in 
P,J(Q) 2 x h N ; for any (t,g) in H A ’(Q) 2 x B A , T A N (f,g) is equal to the function u N , 
where (u N , p N ) is the solution of problem (lll.40) A (ll.27) A in P N (0) 2 x M N . As in the 
continuous case, for f in H A '(0) 2 , T A N f andT AN (f,0) coincide. 

Next, we consider the nonlinear term. Due to ( 1 1 . 1 7 ) A , we need the following 
operator $„ , definedfrom 6°(0) into P N (0) by : for any function f in 0°(O), Sjjf satisfies 
(IV. 9) V ip € P N (0), ( Sj^ f , ip ) A = (f.ip) AN • 

Then, we set for any function w in 0°(O) 2 

(IV. 10) G Ai n( w ) n 2^fri S^(9(3^(WjW))/9Xj - f) and G A N (w) = ( G A N (w) , - g ) . 

This definition is equivalent to 

(iv. 1 1 ) vv n €X n , (g an (w),v n ) a = (L 1 2 = , a(3^(w j w))/ax i ,v N ) AN -(f,v N ) AN . 

Finally, problem ( 1 1 . 1 6) A (ll.9) A has the following equivalent formulation : Finds 
polynomial u N in P^(O) 2 such that 



i 

I 


I 
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OV. 12) a u n + T an G an (u n ) = 0 . 

Problem (11.31 ) A (II. 27) A has the following equivalent formulation : Find a polynomial u N 

in P N (0) 2 such that 

(IV- 1 3) a u N + T a n G a N (u N ) = 0 

Our purpose is now to study these two discrete problems by applying Theorem IV. 1 
respectively to the mappings F N = 1 + T AN G AN defined from Z N = P^(O) 2 into itself and 
F n = 1 + T AN G AN defined from 2 N = P N (0) 2 into itself. 

We know from [M2, Chap. 2, Thm lll,2][BM, Thm IV. 5] that there exist a projection 
operator TT^ from H A (Q) onto P N (0) such that the following estimate is valid for any real 
number s > 1 , 

(IV. M) V<P€H*(0), Ik-T^ipll, A ,o + N ll«P-Tr^tp|| 0 A 0 <cN 1 "*||(p|| $An . 

and a projection operator TT^ 0 from H A0 (Q) onto P^(Q) such that the following estimate 

is valid for any real number s > 1 , 

(IV. 15) | V <p € H^ 0 (O) fl H’(O), 

II »P — N o'P I! 1,A,0 + N II «P — N o'P Ho,A,0 ** C I s ! Il'i’lls.A.fi 
Let us denote by N’ the integral part of (N-D/2. We choose equal to TtjJ, 0 u in the case of 
homogeneous boundary conditions and to TtjJ.u in the case of inhomogeneous ones (this 
definition of seems very complicated, but the fact that u N *u^ belongs to P N .,(0) 2 will 
make the estimates more straightforward, as it will appear later), Due to (IV. 1 4) and 
(IV. 1 5), we note that 

(IV. 16) ||u-u*||, A Ci + N||u-u*|| 0 A ft s:c N 1 ' 9 ||u|| $ A0 . 

The computation of the corresponding constants e N and e N , * N and * N , A n (ti) and A n (tj) 
will be achieved in the next paragraph. 


We begin by stating some results about the linear operators T A N andf AN 
?sit1on IV. 1 : For any t in H A '(0) 2 , the operator T A N satisfies 

ii T <n (‘- v n)a 


(IV. 17) l|T A<N f|| 1A0 <c 


V n €P-( 0) 2 II v N II l.A.P 


(IV- 1 8) lim N _ ll(T A -T AN )f||, A 0 = 0 . 

Moreover, if the solution T A f belongs to H A (Q) 2 for a real number s > 1 , it satisfies the 


estimate 
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(IY.19) ll(T A -T A>N )f|| 1A0 <cN 1 - s ||T A f|| sA0 ; 

If the boundary data flj , J € TU4TL , belong to H a (Tj) 2 for a real number x > 2 and if 
the solution T A (f,g) belongs to H A (G) 2 for a real number s 3* 3 , it satisfies the estimate 
(IV. 20) ll(T A -T AN )(f,g)|| 1 A0 < c ( N * l|T A (f,g)|| sA0 + H Z, j t Z / 4Z I! 0j ll tA rj ) • 


Proof : By Proposition III. I , we obtain at once 

a A,N^ T A,N f ■ v n^ 
lll.A.0 


II t an , H 


A.N' “1.A.0 




sup 


V N « K tA,N 


< c' 


sup 
V N € Pjj(0)‘ 


_0|_,v N ) 

IK 


N ''A 


N'll,A,0 


which is (IV. 1 7). Next, due to the definition of the operators T AN andT AN , the estimates. 
(IV. 1 9) and (IV. 20) have already been stated in (1 1 1.9) and (111.41 ) respectively. Finally, 
(IV. 1 8) holds by classical arguments using (IY.1 9) and the density of £>(0) into H A0 (G). 


In order to estimate the nonlinear term, we need the following lemma. 

Lemma IV.3 : For any real number e > 0, there exists a constant c such that , for any 
y H and <)/ N in P N (0) , the following estimate is satisfied 

(IV. 21) ||(ld-3^)(ip N i|/ N )|| 0 Aft + || ( Id— TC N _ , )(<p N 4 » N ) || 0 A 0 < c N II*Pn IIi,a,o H't'N Hi,a,o • 

Proof : Recalling that N' stands for the integral part of (N-1)/2, we write 

(id-3 N )(ip N »i» N ) = (id-3 N )[(<p N -n N .(p N )(«t< N -n N .i|/ N ) + tt n . <p N (»)/ N -n N . t N ) 

+ ( | Pn - ^N ' ’ 

so that 

||(ld-3(J)(<p N H» N )|| 0 A 0 < ||(ip N -TT N . 4 l N^'t , N — ^ N ' ^Hoao 

+ II ^IIoAO + H^N - Ho,A.O 

+ I|3^[(«P N — TT^. <P n )(4»n — llo,A,0 + II^N fN^N - ^N' *l , N^llo,A,0 

+ II3£ [('Pn-^n 1 1 ) n^n , 'I'n1IIo.a,o 

This implies 

||(ld-3£)(<p N *ii N )|| OAO < c ( ||ip N -TT N . <p N ll|_oo(rt) II'I'n-^n’^n^IIo.a.o 

+ l|TT^. <P N ll L co( 0 ) lk N -Tt^.<|» N || 0 Aft + ||<P N — TT n . <P n || 0 ao l|n N .»|» N || Le »( n ) ) 
^ c ( ll<P N ll L oo (0 , + lln N . ip N || L oo (0 )) l* N - n r*Nll 0A ft 

+ ll<P N - rr N < P N llo.A.ft II ^N‘ Hl°°<0) 

Using the imbedding of H A +t/2 (fi) into L°°(G) for any e > 0 and (IV. M), we obtain for any 

S 3s 1 

(IV. 22) ||(ld-3^)((p N <|/ N )|| OAO < C N 5 [( ||<P N lli +e /2,A,ft + II ' *Pn lll+e/2.A,0 ^ H^N H»,A.O 
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+ H'Pn lls.A.n li ^ N • I1 1 +c/2,A,0 5 ' 

Consequently, we derive the estimate of ||(ld-3jj) (<P N 'V N ) ll 0 A ,n as an aa $y consequence of the 
inverse inequality [CQ1 , Lemmas 2.1 and 2 .4], validfor any integer m and any real number 
r, 0 < m ^ r, 

(IV. 23) Vip N €P N (0), l[tp N ll r ,A,n < c N * II*Pn IL.a.o • 

The term l|(ld-n^_ 1 )(ip N (|» N )|| 0 A0 is estimated exactly in the same way. 

We can now state the following result. 

Proposition IV.2 : For N large enough , the operator DF N (u;j) = 1 + T A N DG A N (u^|) is 
an isomorphism of P^(G) 2 , and x N is bounded by a constant y independent of N. For N 
large enough , the operator DF N (uj$) = 1 + T AN DG AN (u|$)/s an isomorphism of P N (0) 2 , 
and Vn is bounded by a constant ^independent of N. 

Proof : We write DF N (u^J) and DF N (ujJ) in the form 

(IV. 24) DF N (U*)= [1 + T a DG(u)] - (T a -T a n )DG(u) - T A N (DG(u)-DG(u N *)) 

- t a , n (dg-dg an )(u*) , 

and (since the second arguments in G and G A N are constant) 

(IV. 25) DF n (u*)=[1 + T a DG(u)]-(T a -T an )DG(u) -T an (DG(u)-DG(u*)) 

- W d G- dg a,n)( u n*) • 

Since the operator 1 + T A DG(u) is an isomorphism of H A0 (Q) 2 and the operator 
1 + T a DG(u) is an isomorphism of H A (0) 2 , there exists a constant c 0 independent of N such 
that, for any w N in P^(G) 2 , 

(IV. 26) || [1 + T A DG(u)].w^ ||j A ^ > Cq || II f ^ a.o 1 

A 

and, for any w N in P N (0) , 

(IV. 27) II [1 + T a DG(u)].w n II, ao > c 0 || w N ||, Aft . 

It remains to bound the three other terms in (IV, 24) and (IV. 25). Let w N be any polynomial 
in X N . 

1 ) It follows from (IV. 1 8) and from the compactness of the operator DG(u) (see Lemma 
IV. 1) that 

l' m N-.eo IK T A _T A.N^ DG ^ u ^lje.(H A (Q) 2 ,H A (Q) 2 ) ' 

Hence, for N large enough, one has 

(IV. 28) ll(T A -T AN )DG(u).w N || 1AO <(c 0 /4)||w N || 1Art . 

2) It follows from (IV. 17) and from the continuity of the operator DG (see Lemma IV. 1 ) 
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that 

||T aj ,( dg ( u )- d G(u*))||^(h;( 0 )2 H^CO) 2 ) < c || D G(u)-DG(u^) ll i e,(H^(0) 2 ,H^ 1 (Q) 2 ) 

<c'||u-u*|| 1A0 . 

From (IV. 16) together with a density argument, we infer the convergence of ujj to u , 
whence, for N large enough, 

(IV. 29) ||T AN (DG(u)-DG(uJ)).w N || 1A0 <(c 0 /4)||w n || 1ao . 

3) We recall that denotes the orthogonal projection operator from L A (0) onto P N (0) 
and we note that, for any <p N in P N _, (0) , S A ip N is equal to <p N (see (IV. 9)). Thus, by (IV. 3) 
and (IV. 1 0), we know that, for any v N in P^(0) 2 , 

((DG-DG an )(u*).w n ,v n ) 

= L, 2 :) ((a/ax i )(id-n„.,)(w Ni u* + u Ni *w N ),v N ) A 

- ((a/8x.)(3 N -TT N _ ])(w N jll N + U Ni W N^ ’ V N^A,N ’ 

whence, by (IV. 1 7), 

II T a n ( dg - d G a , n ) (u n) w nHi Aft 

< c ( |K 1 d— 3^) ( w Ni uj!j ) llo_A.fi + IKld-n N .,)(w Ni u*) II 0 ,a,o 

+ ll(ld-3 N A )(u Ni *w N ) llo.A.0 + IKld-n A .,)(u N *w N ) || 0 Afl ) . 

From this estimate together with Lemma IV. 3, we derive 

II T A.n( DG-DG A.N^ U nP W N 111, A .ft ^ C 1 II W N Ill.A.n il U N I1 1 ,A,0 ’ 

whence 

(IV. 30) ||T AN (DG-DG AN )(u*).w N ||, Aft ON'" 1 NI, ao ||w n ||, a 0 . 

Finally, we conclude from (IV. 24) to (IV. 30) that, for N large enough, 

V w N € P N ‘(0) 2 , |DF n (u*).w n ||j. A f) ^ K(j/4) I|w n ||, a o , 
and 

V w N € P n (Q) 2 , l|DF N (u*).w N || 1Afl ^(c 0 /4)||w N || 1AO , 

which proves the proposition. 

I ftmrna IV.4 : The constants A N (ig) and Ts^iT)) satisfy 
(IV.31) A n (i))<ct) and A n (-o)^ciq . 

Proof : Let w N be any element in X N . We have 

ll T A,N^ DG A,N^ W N _U N^ll5£,(X N .X N ) 

< ll T A,N ^ DG ^ W N _U N^II't(X N ,X N ) + II T A.n( DG - DG A.n)(w N - u ^) ll;e,(X N ,X N ) • 

Using (IV. 1 7) and the continuity of the operator DG (see Lemma IV. 1) yields that 
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H T A,N (DG( w n — u ^)) ||^(X N x n ) < c II DG ( w N“ u N^I^,(X n ,H^ , ( 0) 2 ) ^ c ll W N“ U N Hi, A, ft 
On the other hand, we know that, for any ^ in P N (0) 2 , 

II T A,N^ DG_DG A,N^ W N _U N) Z Nlll,A,ft 

^ C 2-i il [IK ld -3 N ) [( W Ni -U Np Z N]Ho.A.ft + IK Id - n n - 1 ^ K W Ni -U Ni ) z rJHo,A,ft 
+ || ( Id— 3^ ) [2 Ni ( w N — ) ] II o,a,o + II (Id— TT N _ | ) C z nK w n — u n^ 3 Ho,a,o I ' 
so that, by Lemma IY.3, 

H T A,N (DG-DG A>N )(w N -Ujj) llje,(XN.XN) ^ C N ' ll W N _U N Hl.A.ft 
These two inequalities, together with the definition (IV. I ) of A n (tj) and A n (ti), imply 
(IV. 3 1 ). 

Lemma IV. 5 : The constants e N and e N satisfy 

(IV. 32) e N < c(u) N ,_ * + c(f) u U2x ’' and e N < c(u) N 1_ * + c(f) N U2 “‘ ff + c(fl) N 7/2_t . 

Proof : Using (IV.4) A , we write F N (ujJ) in the form 
= U N + Ta.nGa.n^n*) “ U “ T a G(u) 

= (ujj-u) + (T an -T a )G(u) + T AM (G(U*)-G(u)) + T AN (G AN (U*)-G(u*)) 

which gives 

(IV. 33) e N < llu-u'H, Aft + ll(T A -T AN )G(u)||, Aft + ||T a n (G(u)-G(u n *)) «, ao 

+ ll T A.N^ G ^ U N^ _G A.N^ U N^ I1 1 ,A,0 

It remains to estimate these four terms. 

1 ) Using (IV. 1 6) yields 

(IV. 34) llu-u^H, Afl <cN Us ||u|| s A0 . 

2) It follows from (IV. 19) that 

«< t a- t a.n> g <“> Hi, a .ft <cN'- s l|T A G(u)|| sAn , 
whence, thanks to ( I V.4) A , 

(IV. 35) II(T a -T an )G(u)|| 1ao <cN , - 5 ||u|| $ao . 

3) Due to (IV. 1 7) and to the continuity of G (see Lemma IV. 1 ) , we have 

IIW G (“)- G K>>Hl.A.ft ^ c ll u-u N Hl,A.ft • 

so that 

(IV. 36) l|T AN (G(u)-G(u N ‘))||, Aft <cN , - 5 ||u|| sAft . 

4) From (IV. 17), we derive 
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II t a.n(G( u n)-G a .n( u n))II,ao^ c S 4P.. , 2 

V N € Pn(^) 

Using the definitions (IY.3) and (IV. 1 0) of G and G AN 
Legendre case, 


((G(u^)-G a , n (u^) ,v n ) a 
H v nII).a.o 

, we have for any v N in P^O) 2 , in the 


(G(u*)-G ln (u*),v n ) 

= L ; s1 [(u Ni *u*. av N /ax.) N - (iVu'.avN/axj)] - (f,v N ) + (»,v N ) N , 

and, in the Chebyshev case, 

(G(u*)-G cn (u*),v n ) u 

= Lf =1 [(u N fu N *, (a(uv N )/ax.) u"') uN - (u N *u*,(e(wv N )/ax.) cr’)j 

- (».v N ) u + (».v N ) u N , 

But, since u Ni *ujJ belongs to P N1 (0) 2 , the exact and discrete scalar products coincide. 
Hence, we obtain from (111.37) 


I (G(u^)-G A N (uj)),v N ) A | - I (^.v N ) A - (f,v N ) AN | < cN 




Va,o ,|V n i»iao 


so that 

(IV. 37) ||T a n (G(u^)-G a N (ujJ)) ||, A ft < cN* 2 ll*ll„ A o 

Finally, we derive the desired estimate for e N from inequalities (IV. 33) to (IV. 37). 
The estimate for e N can be obtained in a very similar way, by writing 

F n (u*) = (u*-u) + (T a n -T a )G(u) + T an (G(u*)-G(u)) + T A N (G a n (u*)-G(u‘)) 
and using (IV. 20) instead of (IV. 2 1 ). 


I V • 3 . Existence result and error estimates for the velocitu. 

We can now prove the main results of this section. 

Theorem IV. 2 : Assume that hypothesis (11.20) holds and that there exists a solution 
(u,p) of the Navier -Stokes equations (1. 0(1.4) such that the operator 1 + T A DG(u) is an 
isomorphism of H A0 (Q) 2 ; assume moreover that u belongs to H S (Q) 2 for a real number 
s > I and that the data f belong to H"(G) 2 for a real number o > I . F or N large enough, 
problem (II. 1 6) A (1 1.9) A admits a solution (u N , p N ) in P^(0) 2 x M N . Moreover, it 
satisfies 

( IV. 38) llu-ujl, A 0 < c(u) N'- s + c(f) N U2k -° 
for constants c(u) and c(f) independent of N. 

Theorem IV. 3 : Assume that hypothesis (11.20) holds and that there exists a solution 
(u,p) of the Navier-Stokes equations (I. 0(1.2) such that the operator 1 + T A DG(u) is an 
isomorphism of H A (Q) 2 ; assume moreover that u belongs to H s (0) 2 for a real number 
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s > 3 , that the data f belong to H"(0) 2 for a real number o > 1 and that the boundary 
data gj , J € TL/A7L , belong to H^Tj) 2 for a real number x > 7/2. For N large enough, 
problem (11.31 ) A (II. 27) A admits a solution (u N , p N ) in P N (0) 2 x M N . Moreover, it 
satisfies 

(IV. 39) ||u-u N || 1A0 <c(u)N 1 * s + c(f)N ,+2a - ff + c(g)N 7/2 * T . 

for constants c(u) , c(f) and c(g) independent of N. 

Proof : Using Proposition IV. 2 and Lemmas IV . A and IV. 5, we notice that 2* N A N (2* N e N ) 
and 2i N A N (2* N e N ) are bounded respectively by c e N and c e N ; consequently, the 
assumptions of Theorem IV. 1 are satisfied for N large enough. Hence, there exists a constant 
c > 0 independent of N such that, for each i) < c, there is a unique solution u N of (IV. 1 2) A in 
the ball S N n { w N € P^(O) 2 ; ||w n -Un ||, A0 < i) } (resp. a unique solution u N of (IV. 1 3) A 
in the ball 5 N = { w N e P N (0) 2 ; l|w N -u^||, A 0 ^ x ) }). Next, from (IV. 2), we derive the 
estimate 

K- u n*IIi.a.o < c II f n( u ^Hi.a.o (resp. ||u N -u*||, An ^ c t|F N (u*) || 1>A ,o> - 
which, together with Lemma IV.5, yields (IV. 38) and (IV. 39). 

Next, by Proposition 1 1 1. 1 , there exists a unique p N in M N such that 

^ V N € . blA.N^ V N > Pn) = - a A,N^ U N ' V N^ “ ^®A,N^ U N^ ■ V N^A 

and the pair (u N , p N ) is a solution of the corresponding problem (II. 1 6 ) A (H.9) A or 
(11.31 ) A (ll.27) A . 

Remark IV. 1 : The error bounds we obtain are exactly the same as for the Stokes 
problem ; in particular, the result is still optimal with respect to the regularity of the 
solution (and also of the data f when Chebyshev approximation is used). 


IV.4. Error estimates for the pressure. 

In order to state an error bound for the pressure, we need a lemma. 

Lemma IV. 6 : The approximate velocity u N , as defined in Theorem IV. 2, satisfies 

(IV. 40) sup (G( ^- G a : K ) . v n ^_ ^ c(u) N i-t , c(f) H uu- 

v n €P^(0) 2 H v n Hi.a.o 

The approximate velocity u N , as defined in Theorem IV. 3, satisfies 
(IV.4 1 ) sup (g(u)-G a ,|(u n ),» m )a <c(u ) N '- , c(| ) N I.S»-« ,„(g) N V/S-< 


II V N lll.A.0 


Proof : Let v N be any element in P^(Q) 2 . We consider only the case of homogeneous 
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boundary conditions, since the proof in the general case is strictly the same. We compute 
(G(u)-G an (u n ),v n ) a = (G(u)-G(u n ),v n ) a + (G(u n )-G an (u n ),v n ) a 
Lemma IV. 1 and (IV. 38) give at once 

| (G(u)-G(u n ),v n ) a | < ( c(u) N 1 '* + c(f) n u2 "‘® ) ||v N || ) A 0 . 

From the definitions (IV. 3) and (IY. 10) of Gand G AN , we obtain 
I (G(u N )-G AN (u N ),v N ) A | 

^ c ( =t IK Id - )( u n« u n) Ho,A.n + IK 1 d- TT N _ , ) ( u N i u N ) || 0 , A .o ) H v nII),a.o 

+ I (f.v N ) A - (».v N ) A N | 

Using Lemma IV. 3, we know that, for e > 0, 

||(ld-3 N )(u Ni u N )|| 0 A0 + ||( Id— TT N _ , ) (u Ni u N ) || 0 A 0 

= ||(ld-3 N )(u Nj u N -u N *u^)|| 0 A 0 + ||(ld-n N1 )(u Nj u N -u N *u^; 

= ll(ld-3^)[u Ni (u N -u N *) + u*(u Ni -u Ni *)]|| OAO 
+ |K Id— TT N _ 1 ) [ u Ni (u N — u^J) + u^(u Ni -u N *) 

< C N ( l|u N ||, A n + lludllt A,0 ^ ll U N~ U N Ill.A.O 

< c(u) N *- 1 ( I|u-U N ||, A 0 + ||u— u* II, A 0 ) . 

Then, Theorem IV. 2, (IV. 16) and (111.37) yield for e > 0 

| G(u N )-G A N (u N ),v N ) A | < ( c(u) N' - * + c n 1 * 2 "'* | 

Finally, these two bounds imply (IV. 40). 


"O.A.O 


ll0,A,P 


Hff,A,n ^ H v nIIi,a,o 


Theorem IV. 4 : Assume that hypotheses (11.20) ancf (11.21) hold and that there exists 
a solution (u,p) of the Navier-Stokes equations (l.1)(l.4) such that the operator 
1 + TDG(u) is an isomorphism of H A0 (O) 2 ; assume moreover that it belongs to 
H a (Q) 2 x a real number s > 1 and that the data f belong to H A (0) 2 for a real 

number <x> 1. Then, the solution (u N , p N ) of problem (II. 1 6) A (ll.9) A satisfies 
(IV. 42) ||p-p N || O AO < c(u,p) N 3 -* + c(!) N 3+2R - ff 
for constants c(u,p) and c(f) independent of N. 


Theorem IV.5 : Assume that hypotheses (11.20) and (11.21) hold and that there exists 
a solution (u,p) of the Navier-Stokes equations <1 . 0(1.2) such that the operator 
1 + T DG(u) is an isomorphism of H A (Q) 2 ; assume moreover that it belongs to 
H a (Q) 2 x a real number s > 3, that the data f belong to H A (0 ) 2 for a real 

number a > 1 and that the boundary data g J( J € Z/4Z, belong to H A (rj) 2 for a real 
number x > 7/2. Then, the solution (u N , p N ) of problem (11.3 1 ) A (ll.27) A satisfies 
(IV. 43) ||P-P N lio.A.o < cCu.p) N 3 -* + c(f) N 3+2 “' <r + c(g) N ,,/2 - T 
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for constants c(u,p) , c(f) and c(g) independent of N. 


Proof : Let us introduce the solution (u N ,p N ) in X N x M N of the following problem : 


(IV. 44), 


V V N € ' a A,N^N ■ V N^ + b 1A.l/ V N • P|P + > V N^A “ ® 


^ Pn € ^2A,N * b 2A,N^ N • ® > 

together with the boundary conditions ( 1 1 . 9 ) A (resp. (1 1 . 2 7) A ). Since u N is equal to 
- T a n G(u) (resp. - T an G(u)), we deduce from (IV. 1 9) (resp. (IV. 20)) that 
(IV. 45) Hu-uJI, Afl < c(u) N 1 '* (resp. ||u-u N ||, A 0 < c(u) N Us + c(g) N 7/2 ' T ) ; 

moreover, we obtain from (III. 10) and (111.42) 

(IV. 46) IJp-p N II o.a,o < c(u.p) N 3 ‘ s (resp. ||p-pJ| O AO <c(u,p) N 3 ** + c(g) N 1,/2 ‘ t ) . 
Next, due to ( I V. 44) A and (11.1 7) A or (11.31 ) A , we notice that, for any v N in P^(Q) 2 , 

b 1A, N ( V N -Pn-Pn) = a A,N^N - V n) + <G( U >. V N>A-«A.N< U N - V N> " «W U n) * V N>A • 

so that, from Proposition I II. 1 , we deduce 

||, ^, m 2 a A.N^N-UN. V N) + (G^)-G A . N (u N ),V N )A 

(IV. 47) IIP N -P N ll 0 ,A.n <cN sup jf— ji • 

v N €P*(G) z II v n Im.a.o 

Let v N be any element in P^(Q) 2 . By the uniform continuity of a A N , we have 

a A,N^ u N“^N ■ V N^ ^ C K^N Hi, A, ft H v N Hi, A, ft 

< C ( I|u-U N ||, A 0 + I|u-U N ||, A n ) ||v N ||, A 0 , 

so that one can bound this term from (IV. 38) or (IV. 39) and (IV. 45). Using this estimate 
and Lemma IV. 6 in (IV. 47) yields 


,3+2k-it 


(IV. 48) lfP N -P N II 0 , A ,0 < c ( u ) N + °( f ) N 

(resp. II P N -P N llo.A.n < °( u ) n3 ’ $ + c < f ) N 
which, together with (IV. 46), gives (IV. 42) and (IV. 43). 


3+2k " + c(g)N 1,/2 - T ) 


That ends the theoretical results which can be proven for both the Legendre and 
Chebyshev approximations of the Navier-Stokes equations. It remains to apply this method 
to real problems, as will be done in the next section. 
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In this section, we describe the resolution algorithm we use for numerical 
applications. It has been proposed first by Y. MORCHOISNE [Mo] and is aimed at solving the 
time-dependent Navier-Stokes equations 

(Y.1) 0u/dt - v Au + grad p + (u.V)u = f inQx(OJ) ,T>0 , 

div u = 0 in 0 x (0,T) , 

with initial condition u(0) = u° in 0. But it can also be used to compute stationary cases as 
it will be shown in the following. As far as tkne-dependent problems are concerned, time 
discretization is achieved with the help of a finite difference scheme. While the convection 
term is handled explicitly by an Adams-Bashforth approximation, the diffusive term is 
implicitly treated in order to ensure stability. 

Let us introduce a fixed time step St > 0. At each time (n+ 1 ) St, n > 0, we compute an 
approximation uf^ 1 in X N of the velocity u((n+ 1 ) St). Furthermore, in order to make the 
numerical computation easier, we first compute a scalar quantity q^*' in P N (0), that we 
call the pseudo-pressure, such that grad q£ +1 is an approximation of the pressure gradient 
(grad p)((n+1) St). When the convergence is reached, the discrete pressure p N is then 
obtained by a post-treatment which is performed by solving a Poisson problem. 

Numerical applications (see § Y.4) have been made with a Chebyshev spectral 
discretization. Thanks to this choice, we can employ the Fast Fourier Transform (FFT) in 
the computation of the derivatives (see [CLW][CT][02]). 




We consider the Navier-Stokes equations (Y. 1 )(l.2) with null right-hand side f. For a 
given function g satisfying the assumptions (1.15), (1.16) and (11.26), we introduce the 
subspace X N (g) of all polynomials in X N satisfying the boundary conditions (ll.27) A . 

Let (u^ , q,J) be any initial quantities in X N (g) x P N (0). We assume that (uj] , qJJ) is 
known in X N (g) x P N (0), and we seek (uJJ +1 , q£ +1 ) in X N (g) x P N (0) such that 

u n+1 -u n 

[L N N v AuJ]* + grad q^’ + (uJ]*.V )u^*](x) = 0 , x € H 0 , 

ot 

(V.2) 1 1 

(div u|J + , div ) A N = inf (div W N , div w N ) A N . 

w N € X N (g) 

In the equations (V.2), we use the following notation : for any integer n > 1 , 
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u„* = (3/2) uJJ- (1/2) u ^' 1 . 

Moreover, the operator L is a finite difference approximation of Id - i)A, where i) is a 
positive parameter. More precisely, we set L = L,L 2 , where 
(V.3) L.rld— Tjd, , i = 1 or 2, 

and dj is the second-order finite difference operator : for example, if i is equal to I , for any 
function w in ($°(Q), we define for any nodexj A = (C A , Cjf) in fl 0 



2 




X 


r • *?) - C k A ) w(C A , C k A ) - w(cf_, . C k A ) , 

t A -t A r A -r A ' 

S S S-i 

The parameter i) verifies 
-g = § v 6t + y V 8t 2 , 

where g and y are two nonnegative constants, and V is an estimate of the velocity norm. Note 
that, if L is chosen equal to the identity (i.e. , i) = 0) , we have an explicit Adams-Bashforth 
scheme. In fact, however, we choose & and y large enough to ensure a good stability of the 
scheme. Indeed, this scheme, when applied to the one-dimensional Burgers' equation, has 
been analysed in the periodic case; it has been proven [Me, Chap. 1 ] that for i) large enough, 
it is unconditionally stable and has a precision upper-bounded by c(a) (v 6t + 6t 2 + N' ff ) 
for all real numbers a > 0. 


V.2. Velocitu and pseudo-pressure computation. 

Problem (Y.2) is solved in two steps. First, we compute a predictor uJJ*p of the 
velocity in x£(g). Then a corrector (v N , q N ) in P^(0 ) 2 x P N (0) is computed, so that the 
pair (ujj* 1 , qJJ* 1 ) defined by 
(V.4) I <’=<;.», . 

verifies the equations (V.2). 


(i) Velocitu predictor computation. 

We first solve the following problem : Find a polynomial u^p in XjJ(g) such that 
u-i u » 

(V.5) [L -^ p N - v AuJJ* + grad q£ + (u^*.V)u^](x) =0 , x € 2 ( fl 0 . 

ot 

This problem can be handled with standard linear system algorithms. Indeed we can 
associate with each operator L, , i = 1 or 2, the operator L, defined from G°(0) into P N (0) 
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bg : for any function w in 6 °( 0 ), Lw belongs to P N (0) and satisfies 


(V. 6 ) 


LjW(x) = LjW(x) ,x« ojflfi , 

L i w(x)nw(x) ,X€ fiaO . 

From [V, Thm VI. 3], we deduce that the operators L. , i = 1 and 2, are (easily) invertible in 
P n (Q) for any 13 > 0. We set L = L,L 2 . 

Let s be in X N , the solution w in X N of the problem Lw = s is obtained by solving 
successively the two following problems : find w, in X N such that CjW, = s and find w 2 in X N 
such that C 2 w 2 = w, ; thus we have w = w 2 . 


(ii) Velocitu corrector and pseudo- pressure computation. 

Thanks to (V.4), the pair (uJJ +1 , q^* 1 ) is the solution of (Y. 2 ) if and only if the pair 
(v N , q N ) of P N '( 0) 2 x P N (0) satisfies 

(L + grad q N )(x) = 0 , x € z£ fl 0 , 

(V.7) 

(div u N n+ p + div v N , div u; + p + div v N ) A N = inf (div w N , div w N ) . 

w N e X*(g) 


In order to solve the problem (V.7), we introduce an operator from P N (0) into 
P N (0) which connects the pseudo-pressure q to div v. We first define the operator grid in 
the following way : for any r N in P N (0), grid r N belongs to P N (0 ) 2 and satisfies 
(V. 8 ) (grid r N )(x) = (grad r N )(x) , x € ZjJ D 0 , 

( grid r N )(x) = 0 , x € ZjJ fl 60 . 

Then, we set, for any r N in P N (0), 

(V.9) ^r N = - 8 t div (L _1 grid r N ) . 

Thus, we can consider the two following problems : Find q N in P N (0) such that 


(V. 1 0) (<J&q N + div u^ + p , ^q N + div u{J*p) A _ N 

= in i / \ ( ^ r N + d1v u n. + p ’ ^ r N + div 
r N € P„(Q) 

and, secondly : Find v N in P^(Q) 2 such that 

(V. 1 1 ) (L + gradq N )(x) = 0 , x € Z* fl 0 . 

01 

In the computations, we handle the resolution of (Y.7) by solving the system (V. 1 0) ( V. 1 1 ). 
Hence, we are going to prove that the equations (V. 1 0) and (V. 1 I ) are equivalent to (Y.7). 


First, we need some results in order to prove that the minimization problem (V. 1 0) 
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is equivalent to the minimization problem in (V.7). 

Lemma V. 1 : If the parameter T) is small enough, the kernel of the operator is equal 
to Z 1A N and its range is equal to M 2A N , 

Proof : Clearly, Z )A N is contained in the kernel of and the range of <$$- is contained in 
the range of the divergence operator, hence in M 2A N . Consequently, it suffices to prove that 
the kernel of S& is contained in Z )A N , since that this would imply that it is equal to it and of 
dimension 8, and that the image of ^-andM 2AN have the same codimension 8 in P^fi) 2 . 

Thus, let q N be any polynomial in the kernel of S&. That implies that C" 1 grad q N is 
divergence-free in 0; since it belongs to P^CO) 2 , there exists a unique polynomial >|> N in 
P N+ ,(0) fl H 2 0 (O) such that 
(Y.12) L* 1 grad q N = curl <j/ N in 0 

Writing the expansion of i|» N in the form 

4 > n (x) = ( 1 -x 2 ) Ln-, 'kn(y) )’(x) , 

we obtain (5(>k N ( 1 -x 2 ) K )/9x)(x) = t n (y) n(n+2oc + 1 ) J*(x) ( 1 -x 2 ) K , so that the 
degree of <j/ N with respect to y is < N. Using a similar argument for the variable x, we 
deduce that >|> N belongs in fact to P N (0). Next, we compute 

( grid q N , curl (^ N ( 1 -x 2 ) K ( 1 -y 2 )*) ( I -x 2 )‘“( I -y 2 )"“) A N 

= (grad q,g , curl( 4 r N ( 1 -x 2 )“( 1 -y 2 )*) ( 1 -x 2 )'*( I -y 2 )*") A N 
= j 0 grad qN . curl («|/ N ( 1 -x 2 ) K ( 1 -y 2 )*) dx = 0 . 

That implies by (V. 1 2) 

(LCcurl 4 » n ) , curl (* N ( I -x 2 ) K ( 1 -y 2 )“) ( 1 -x 2 )'"( 1 -y 2 )' K ) AN = 0 , 

or equivalently 

(LCcurl y m ) , curl (* M ( 1 -x 2 )“( 1 -y 2 )“) ( I -x 2 )'“( 1 -y 2 r*) AN = 0 . 

1 ) In the case i) = 0, we have proven that 

((curl »n N ) , curl (<j, N ( 1 -x 2 )“( I -y 2 )*) ( 1 -x 2 )-“( I -y 2 )* K ) A , N = 0 . 

From the ellipticity of this form on Pn( 0), we deduce at once that qi N is equal to 0, hence 
that grid q N is equal to 0, and q N belongs to Z )A N . 

2) In the case t> > 0, denoting by c(N) the norm of the operators 9, and 0 2 on the space 
P N (0) provided with the norm ||.||, A and writing L - id = - r) 9, -i} 9 2 + -g 2 9,9 2 , we 
obtain 

0>c ||y n II 2 A - ti [2 c(N) + t) c(N) 2 ] ||.|< N || 2 a0 , 
whence the result for ig small enough. 


! 
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From the two results of this lemma, we derive respectively the two following 
propositions. 

Proposition V. 1 : If the parameter 13 is small enough, the system (V.IOHV.ll) is 
equivalent to problem (V.7). 

Proposition V.2 : If the parameter i) is small enough, the set of values {(grad q N )(x) , 
x € 3 N A fl 0 } , where q N is a solution of problem (V. 1 0 ) , is uniquely defined and the 
solution v N of problem (V. 1 1 ) is uniquely defined. 

In both Legendre and Chebyshev cases, the minimization problem (V. 1 0) can be solved 
thanks to the Axelsson's minimization algorithm, which was aimed to problems associated 
with symmetrical nonnegative operators or with operators the symmetrical part of which 
is positive definite [Ax][J][Me]. In our case, even if the operator does not satisfy these 
assumptions, the algorithm turns out to be efficient when appropriate re-initializations 
are used [Me, Chap. 2, § VI 1.3 and Chap. 4, § IV. I ]. 

Remark V.l : Note that it is rather standard to set up problems concerning pseudo- 
pressure, as in (V. 10), by eliminating the velocity of the continuity equation. The basic idea 
of this procedure relies upon the Uzawa's algorithm [61], since the pseudo-pressure plays 
the role of the Lagrange multiplier. 


V.3. Pressure post-treatment. 

Once we have reached the stationary state, i.e. , the time (n + 1 ) 8 t when the velocity 
becomes independent of the integer n, we can compute the pressure. 

We set 

u n+, -u n 

(Y.13) S "* 1 = L N N vAu^tuJ'.V)^* , 


and we seek the pressure p M in P N (0) as the solution of 
(V.l 4) Ap N (x) = - (div S n 4 , )(x) ,x(Z N A flfi . 

(0p M /0n)(x) = -(S nt 1 .n)(x) , x € Z* f) 90 . 

Here, the vector n denotes the unit outward normal vector to 0 on 90. 


Problem (Y.14) is solved through a finite difference preconditioning method which 
involves the operator L defined in (V.3). Thus, p N is computed as the limit of the following 
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sequence ( p Nk ) k2;0 . We set p N0 = q^’ ; then, assuming that p Nk is known, we compute 
p Nk+1 as the solution in P N (0) of 

(V.15) (L(p Nkt1 -p Nk ))(x) = X (Ap Nk + div S n+, )(x) , x € D 0 , 

^ 9 ^PN,k+i - PNp /an |DF^ x ^ = “ y (Op Nk /0n + sn+ ' n )( x ) , x <E 3* f) 00 , 

where X and y are suitable parameters and the operator a/9n| DF is defined as follows : for 
any polynomial r N in P N (0), we set for each point x of fl 90 , 

1 ) if x is not a corner, assuming for instance that x =x ok = (Co * ) belongs to r, , 

r N • ^k ^ “ r N ^ 1 > ^k ) 


(0r N /en |oF )(Xo k ) = 


-pTTF* 

‘♦o ^ 1 


2) if x is a corner, assuming for instance that x is equal to a, = (C© , O' 

«Va V )(V . . 

the operator 0/dn| DF is defined similarly on the three other edges l"n , r H | and T| V of the 
square. 

The parameters X and y are chosen experimentally, in order to ensure the 
convergence of the sequence ( p Nk ) k?0 . In our computation, they are respectively equal to 
0.01 and 0.75. 


The finite difference preconditioning method is well-known for spectral computations. 
It allows one to avoid direct inversion of spectral operators (e.q. the operator A in our 
case), which is expensive because the corresponding matrices are full (see [02], for 
instance). 

Finally, let us remark that we did not look for the pressure in the space M N as is 
suggested in the theory of § 1 1 1.3. This approach is now under consideration. 

V.4. Numerical results. 

The numerical experiments were performed in domains of [R 3 with curved geometries, 
in which we generalized the previous algorithm. Indeed, we can use spectral techniques in a 
curved connected open set 0 c IR 3 , if there exists a one-to-one function & which is 
sufficiently smooth and maps the reference cube 0 = ]- 1 , 1 [ 3 onto 0. Thanks to the function 
& , a problem initially set in 0 is brought back to the cube 0. 

The set of spurious modes for pressure can be identified for three-dimensional 
problems [BMM , Lemme V. I ] [ B CM , Remark IV. 3], and we can obtain a well-posed problem 
for the Navier-Stokes equations in a cube as has been done in the two-dimensional case. 
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Moreover, we refer to [Me, Chap. 2 and 3] for details about discretization of Navier-Stokes 
equations in curved geometries. 

Let x = (x, , x 2 , x 3 ) and x = (x, , x 2 , x 3 ) denote the generic points of 0 and 0, 
respectively. In our numerical applications, we have considered the set 0 to be defined as a 
curved hexahedron of IR 3 , two opposite sides of which are plane and parallel. Without 
restriction, we assume that these two sides are parallel to the plane x 2 = 0. The nozzle 0 is 
then defined from its boundary : let : A = [ — 1 , 1 ] 2 — ► IR 3 , 1 < i < 4 , be the 
parametrizations of the four other sides (see Figure V. 1 ). 



The function S' can be defined as follows : for x = (x, , x 2 , x 3 ), we set 
^(x) = (( 1 +x 3 )/2) £j(x, , x 2 ) + (( 1 -x 3 )/2) f,(x, , x 2 ) 

+ ((1 +x,)/2) [&(x 2 , x 3 ) - (( 1 +x 3 )/2) £ 4 (x 2 , +1) - (( J -x 3 )/2) f 4 (x 2 , -1)J 
+ (( 1 -x,)/2) [f 2 (x 2 , x 3 ) - (( 1 +x 3 )/2) f 2 (x 2 ,+!)-((! -x 3 )/2) f 2 (x 2 , - 1 )] . 
Clearly, we have S'(dQ) = 90. Moreover, we assume that the function & is one-to-one 
and that S'i 0) = 0. This last property can be deduced from hypotheses of smoothness 
concerning the function S' (see [Me, Chap. 3, Th. V. 1 ]). 

In the example below, the open set 0 is a nozzle, the cross section of which in any 
plane x, = constant is a rectangle. The function S' is simply defined by 
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^(x) = 


x, = a x, + 0 


X 2 = X 2 


X 3 ~ X 3 

where £ is a fourth-degree polynomial with maxima at x, = ± 1 and a minimum at a point 
x, = x° of ]- 1 ,1 [. The boundary conditions g are given by 


, x, = ±1 


3 (l-Xj) x 3 £'(x,) / 4<x£(x,) 



Note that we can extend the boundary data g into a divergence-free function u b on 0, by 
setting 

u 

The viscosity coefficient v is equal to 1 O' 2 . The discretization is performed in the space of 
polynomials of degree < N, with N = 16. Thus the mesh is made up of 17 3 collocation nodes 
associated with the Chebyshev polynomials. The time step 8t is equal to 1 0' 3 . 


, 3 ( 1 -x 3 ) / 4 £(x,) 
o F(x) = | 0 ! r 1 

3 (1-Xj) x 3 £'(x,) / 4oc^(x,) 


g o ^(x) = 


3 O-xf) /4^(x,) 
0 

3 ( 1 -x|) /4^(x t ) 


Figure V.2 shows the mesh. In Figure Y.3, the velocity iso-norms and the 
iso-pressure lines are presented in the plane x 2 = 0, at the time 1 0 6t (when the stationary 
state is already reached). 

Figure V.4 shows the iso-pseudo-pressure lines. The spurious mode T N (x,) (the 
extrema of which coincide with the vertical lines of the mesh) appears clearly and totally 
hides the pressure behavior. 

Figure Y.5 shows the convergence of the algorithm (V.13) for the pressure 
post-treatment. Due to Neumann boundary conditions, the convergence is rather slow, so 
that the technique must be improved in order to obtain an efficient pressure solver. 

Nevertheless, we obtain good results concerning the velocity. Note that we have 

((div u N n+1 ) o & , (div ujj* 1 ) » ^F)^ - 6.4 . 10' 4 . 

This quantity can of course be reduced by increasing the number of collocation points. We 
refer to [Me, Chap. 4, § IV. 3] for another way to reduce ((div uJJ +1 ) o ^ , (div u{] +1 )o^F) C N . 



• 5 - 


Resldual at the Internal nodes 

Boundary conditions at the internal nodes of the faces 


Number of Iterations 
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